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ABSTRACT 

A  modern  computer  based  automatic  data  acquisition/con- 
trol system  was  installed  at  the  Department  of  Aeronautics' 
Combustion  Laboratory.   This  system  utilizes  an  HP-85  desktop 
computer  as  system  controller  for  the  HP-3054A  data  acquisi- 
tion system.   These  major  control  applications  include: 
(1)  particle  size  determination  by  measurement  of  scattered 
laser  light,   (2)  a  vitiated  air  heater  system,  and   (3)  a 
solid  fuel  ramjet  test  unit.   Details  involving  specific  data 
acquisition/control  techniques  for  these  applications  are 
discussed. 
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I.   INTRODUCTION 

This  thesis  project  was  undertaken  to  develop  and  install 
a  modern,  computer  based  experiment  control/data  acquisition 
system  at  the  Naval  Postgraduate  School  Aeronautics  Depart- 
ment, Combustion  Laboratory.   Previously,  all  laboratory 
experiments  were  manually  controlled  and  all  test  data  were 
reduced  manually  by  hand  with  no  real  time  data  processing 
possible.   Precise  timing  of  test  equipment  and  control 
features  were  difficult  or  impossible. 

The  introduction  of  the  HP-3054A  Autom.atic  Data  Acquisi- 
tion system  and  the  integration  of  all  associated  test  equip- 
ment to  a  desktop  computer  have  revolutionized  current 
operations  at  the  Combustion  Laboratory. 

Large  quantities  of  high  quality  data  can  be  acquired  at 
high  data  rates  and  stored  for  future  use  or  processed  real 
time  in  closed  loop  feedback  control  circuits  for  precise 
process  control  and  display. 

Other  electromechanical  equipment  such  as  Visicorders 
provide  analog  real  time  data  and  run  in  parallel  with  digi- 
tal equipment  but  can  be  remotely  controlled  for  precise 
timing  of  experiments. 

Interactive  programs  on  the  HP-85  computer  allow  the  ex- 
perimenter to  do  preliminary  calculations,  initial  experiment 
set-up,  and  to  perform  accurate  process  control  of  system 
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functions.   On  line  data  reduction  and  display  of  results  is 
performed  by  the  system  CRT  display,  the  7225B  Plotter,  or 
the  system  printer. 

Four  basic  experiments/experimental  apparatuses  are  cur- 
rently automated  using  the  automatic  data  acquisition  system; 
a  light  scattering  experiment,  a  vitiated  air  heater,  a  solid 
fuel  ramjet,  and  a  turbojet  combustor  test  rig.   The  first 
three  applications  are  discussed  in  subsequent  sections  along 
with  the  measurement  techniques  used  to  implement  the  process 
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II.   SYSTEM  AND  COMPONENT  DESCRIPTIONS 

To  produce  accurate  test  results  involving  high  speed 
machinery  and  processes  and  simultaneously  provide  precise 
timing  and  control  signals,  a  modern  data  acquisition  system 
has  been  installed  at  the  Combustion  Laboratory  at  the  Naval 
Postgraduate  School. 

Figure  1  shows  a  system  block  diagram  of  the  Hewlett- 
Packard  3054A  Digital  Data  Acquisition  System.   The  heart  of 
this  system  is  the  HP-35  microcomputer  functioning  as  system 
controller,  although  any  computer  compatible  with  the  Hewlett- 
Packard  Interface  Bus  (HPIB)  can  be  substituted  or  connected 
to  the  system  bus. 

A  description  of  each  equipment  connected  to  the  bus  is 
described  below. 

A.   3054A  DATA  ACQUISITION  AND  CONTROL  SYSTEM 

The  3054A  is  a  computer-based  automatic  data  acquisition 
and  control  system.   The  system  is  complete  with  a  system 
mainframe,  the  HP-3497A,  and  two  voltmeters  (the  HP-3456A 
and  HP-3437A) ,  that  are  interfaced  via  the  82937A  HP-IB  I/O 
card  to  the  powerful  HP  desk  top  computer,  the  8  5A. 

The  3054A  System  documentation  has  been  designed  to  co- 
ordinate efficient  use  of  each  of  the  system  instruments, 
and  it  includes  a  specialized  software  support  package  that 
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Figure  1.   Data  Acquisition  System 


is  designed  to  use  the  System  instruments  to  achieve  specific 
functions  and  operations. 

B.   HP-3456A  DIGITAL  VOLTMETER 

The  3456A  Digital  Voltmeter  is  a  3%  to  Sh   digit,  inte- 
grating voltmeter  with  relatively  high  speed,  10  parts  per 
million  (PPM)  basic  accuracy  and  100  nanovolt  sensitivity. 
Its  DC  and  True  RMS  AC  voltage  maximum  of  1000  volts  and  re- 
sistance measurement  capability  reduce  the  amount  of  signal 
conditioning  necessary.   It  can  detect  100  nanovolt  changes 
in  100  millivolt  signals  at  speeds  of  48  readings  per  second. 
This  capability  is  required  for  measuring  thermocouples  with 
the  3054A  system  to  better  than  Q.1°C   resolution.   By  se- 
lecting Ah   digits  of  resolution,  reading  rate  is  increased  to 
330  readings  per  second.   The  3456A  has  selectable  integra- 
tion times  from  0.01  to  100  power  line  cycles  and  provides 
input  guarding  to  give  140  db  of  common  mode  rejection,  pro- 
viding for  accuracy  in  the  presence  of  noise.   Two  other 
noise  reduction  features  are  provided  by  digital  averaging 
and  an  analog  input  filter. 

Up  to  350  readings  can  be  stored  in  the  3456A  Read/Write 
memory  in  addition  to  providing  a  programmed  time  delay  be- 
tween readings.   A  "voltmeter  complete"  signal  provides  the 
necessary  means  to  synchronize  the  operation  of  the  3456A 
with  the  3497A  to  allow  analog  scanning  rates  up  to  300 
channels  per  second. 
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C.  3437A  DIGITAL  VOLTMETER 

The  3437A  is  a  high  speed,  3k   digit,  DC  voltmeter  which 
provides  precisely  timed  sample  and  hold  readings.   It  can 
perform  a  variety  of  tasks  from  scanning  DC  inputs  to  AC 
waveform  analysis.   Repetitive  signals  with  frequency  com- 
ponents up  to  1  MHz  and  low  frequency  transients  longer  than 
500  ysec  can  be  rapidly  digitized  and  analyzed.   Parameters 
such  as  RMS  value,  residual  DC,  harmonic  content,  and  peak 
values  can  be  readily  and  automatically  determined. 

Using  the  HP-85  controller  and  a  "Fast  Handshake"  trans- 
fer mode  of  operation,  reading  rates  of  up  to  4000  readings 
per  second  are  possible  in  the  packed  binary  mode  of  operation 

D.  3497A  DATA  ACQUISITION/CONTROL  UNIT 

The  3497  data  acquisition/control  unit  combines  precision 
measurement  capability  with  control  functions.   It  can  provide 
precision  measurements  of  strain  gauge  outputs,  thermocouples, 
pressure  transducers,  and  other  sensors,  and  can  also  provide 
digital  interfacing  and  control.   A  digital  clock  is  built 
into  the  3497A  mainframe  for  real  time  control  applications. 

Five  slots  are  available  for  plug-in  I/O  card  options  on 
the  3497A  mainframe.   Additionally,  the  3497A  capabilities 
may  be  expanded  through  the  addition  of  one  or  more  3498A  Ex- 
tenders.  Up  to  1000  analog  channels  and  1300  digital  channels 
can  be  obtained  with  the  addition  of  14  3498A  Extenders.   Each 
extender  has  10  slots  for  plug-in  I/O  card  options. 
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The  3497A  mainframe  also  provides  several  key  features 
that  make  it  very  suitable  for  use  in  computer  controlled 
data  acquisition  systems.   The  following  TTL  compatible  input 
and  output  terminals  are  standard  on  the  3497A:   external 
trigger  input,  external  increment  input,  break  before  make 
(BBM)  sync  pulse  (input /output) ,  voltmeter  complete  output 
signal,  channel  closed  output  signal,  and  a  timer  signal. 
Additionally,  the  3497A  has  multiple  programming  commands 
providing  versatility  in  process  control  applications. 

E.  34  98A  EXTENDER  UNIT 

The  capacity  of  3497A  mainframe  may  be  expanded  through 
the  addition  of  from  one  to  fourteen  3498A  extender  units, 
each  having  10  slots  in  its  frame  for  I/O  modules  of  various 
types  to  meet  specific  application  requirements. 

Special  attention  to  the  3498A  extender  operating  manual 
is  required  when  adding  additional  I/O  cards.   Slot  numbers 
5  through  9  do  not  exist  for  digital  cards. 

F.  HP-82901M  FLEXIBLE  DISC  DRIVE 

The  mass  storage  requirements  for  the  system  are  handled 
by  the  8290 IM  dual-drive  disc  system  which  supports  two  5^ 
inch  flexible  double-sided,  double  density  discs,  providing 
a  total  of  approximately  550K  bytes.   Drive  #0  is  established 
as  the  default  mass  storage  medium.   Over  220K  bytes  of  mass 
storage  are  also  available  in  the  HP-85  internal  casette  tape 
drive. 
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G.   HP-7225B  PLOTTER 

The  HP-7225B  plotter  contains  a  17601A  personality  module 
and  is  therefore  a  microprocessor  based  HP-IB  plotter  that 
produces  high  quality  graphic  plots  on  any  size  chart  up  to 
210  X  297  mm. 

Pen  movement  is  as  small  as  .032  mm  for  high  resolution 
plotting.   The  plotter  instruction  set  includes  39  different 
instructions  to  provide  such  capabilities  as  point  digitizing, 
labeling,  character  sizing,  scaling,  and  window  plotting.   The 
7225  is  interfaced  through  the  HP-IB  providing  extreme 
versatility. 

Pen  velocity  is  programmable  from  10  mm/sec  to  250  mm/sec. 
Seven  different  dashed-line  formats,  symbol  mode  plotting  and 
user  defined  characters  aid  in  trace  identification. 

H.   HP-2631B  LINE  PRINTER 

The  2631B  is  a  180  character  per  second  printer  that  pro- 
duces high  resolution  dot  matrix  printing.   Characters  are 
printed  from  one  of  two  128  ASCII  character  sets.   A  16  chan- 
nel fixed  or  programmable  vertical  forms  control  (VFC)  is 
available  for  any  desired  page  and  text  length  with  several 
print  pitches  to  select  from.   An  underline  mode  is  also 
available. 

The  2631B  uses  a  bidirectional  printing  head  and,  addi- 
tionally, monitors  incoming  data  to  determine  optimum  print 
direction  for  the  next  line.   Other  features  include  blank 
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space  detection  and  suppression,  eight  print  densities  with 
underlining,  automatic  vertical  and  horizontal  tab  settings 
and  margins.  It  also  includes  an  end  of  line  (EOL)  wrap- 
around feature.  When  the  printer  receives  data  which  would 
cause  printing  to  exceed  the  right  margin,  it  automatically 
moves  the  data  (which  would  ordinarily  be  lost)  to  the  next 
line. 

I.   HP-85A  COMPUTER 

The  heart  of  the  Data  Acquisition  system  is  the  85A  desk- 
top computer  which  functions  as  system  controller.   The  com- 
puter is  interfaced  with  system  components  via  the  82937 
HP-IB  I/O  card.   The  85A  uses  enhanced  BASIC  programming 
language  and  supports  the  PLOTTER/PRINTER  ROM,  INPUT/OUTPUT 
ROM,  MATRIX  ROM,  and  the  MASS  STORAGE  ROM,  providing  an  ex- 
tremely versatile  instruction  set.   The  85A  also  features  a 
32x16  CRT  display  with  full  graphics  capability,  an  internal 
32  column  thermal  printer,  live  keyboard,  buffered  I/O,  pri- 
ority intercept,  high  speed  tape  cartridge,  autostart,  and 
error  trapping. 

It  is  also  a  versitile  calculator  and  can  be  used  with 
virtually  no  previous  HP-85  experience.   The  user  can  halt  an 
executing  program  using  the  pause  key,  perform  a  series  of 
calculations  on  the  keyboard,  print  out  the  answer,  and  then 
resume  program  execution  by  pressing  the  CONTINUE  key. 

Reference  1  describes  the  HP-85  system  characteristics  and 
the  extensive  instruction  set. 
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J.   HP-6942A  MULT I PROGRAMMER 

For  high  speed  data  acquisition  applications  the  HP-3054A 
is  augmented  with  the  6942A  multiprogrammer .   The  6942A  is 
capable  of  functioning  independently  from  the  30  54A  as  a 
master  control  unit  for  bidirectional  communication,  data 
processing,  and  data  transfer  between  an  HPIB  controller  (the 
HP-85  computer)  and  a  multiprogrammer  system.   It  can  be  used 
in  a  single-unit  system  employing  from  one  to  sixteen  plug-in 
I/O  cards,  or  in  a  multi-unit  system  consisting  of  one  6942A 
master  unit  and  up  to  seven  6943A  extender  units.   Each  ex- 
tender unit  can  accommodate  up  to  sixteen  input/output  (I/O) 
cards,  allowing  a  multiprogrammer  system  to  be  expanded  to 
128  I/O  cards. 

The  6942A  employs  32  different  instructions  in  its  in- 
struction set,  permitting  a  wide  variety  of  I/O  card  func- 
tions under  program  control  of  an  HPIB  controller.   The 
instruction  set  includes  system  control,  output,  input,  card 
control,  and  system  timing  functions.   System  control  in- 
structions establish  the  basic  operating  modes  of  the  multi- 
programmer  system.   Output  instructions  are  used  to  send  data 
to  output  cards  while  input  instructions  are  used  to  obtain 
data  from  input  cards.   The  card  control  instructions  allow 
access  to  various  circuits  on  the  I/O  cards  to  enable  close 
control  and  monitoring  of  their  various  operations.   The  sys- 
tem timing  instructions  permit  the  sequencing  of  events  and 
the  measuring  of  elapsed  time,   A  detailed  description  of  the 
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instruction  set  is  described  in  the  6942A  User's  Guide  [Ref. 
2]  as  well  as  programming  details  for  the  HP-85  controller. 

In  addition  to  the  16  I/O  cards,  an  HP-IB  interface  board 
and  a  transmission  system  board  are  installed  in  the  rear  of 
the  6942A  unit.   The  6942A's  internal  mainframe  circuitry  in- 
cludes a  CPU/ROM  board  and  a  control/RAM  board.   All  power 
supplies  for  the  I/O  card  slots  and  mainframe  circuitry  are 
built  into  the  6942A  unit.   The  CPU/ROM  board  includes  a  16 
bit  microprocessor  chip  which  decodes  and  executes  all  in- 
structions, and  controls  all  data  transfers.   The  12K  x 16  bit 
ROM  (Read  Only  Memory)  contains  the  microprocessor  controlled 
programs  (firmware)  required  to  process  all  instructions. 
The  ROM  also  contains  diagnostic  programs  for  built-in  self 
test  and  signature  analysis/troubleshooting  capabilities. 
The  control/RAM  board  contains  a  2K  x  16  bit  FJ\M.    (Random  Ac- 
cess Memory)  which  is  used  for  the  temporary  storage  of  in- 
structions and  data. 

Backplane  control  circuits  on  the  control/RAM  board  de- 
code I/O  card  addresses  and  control  I/O  card  functions.   The 
HP-IB  interface  board  provides  the  bidirectional  communica- 
tions interface  between  the  6942A  and  an  HP-IB  controller. 

Three  I/O  card  types  are  presently  installed  in  the  6942A 
multiprogrammer :   a  69736A  Timer/Pacer  card,  two  69751A  high 
speed  analog  to  digital  converter  cards,  and  two  6  9790B  4K x 
16  bit  Random  Access  Memory  cards  . 
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1 .   69736A  Timer/Pacer  Card 

At  power-on  the  timer/pacer  "wakes  up"  in  the  one 
shot  mode  which  produces  a  single  output  pulse  with  a 
crystal-controlled  duration  programmable  from  one  micro- 
second to  about  eighteen  hours  (65535  seconds) .   Another 
mode  (the  continuous  mode)  can  also  be  programmed  in  which 
the  output  goes  high  for  the  programmed  time  interval,  then 
low  for  the  same  interval,  and  continuously  repeats  this 
cycle  until  the  card  is  reprogrammed  or  disabled  externally. 
The  card's  single  output  pulse  can  be  used  as  a  programmable 
delay  or  as  an  enable  signal  for  a  frequency  counter  card 
(69775A)  in  frequency  measurement  applications.   The  square 
wave  output  generated  in  the  continuous  mode  can  be  used  to 
pace  A/D  or  D/A  converters  as  they  analyze  or  generate  wave- 
forms or  to  drive  frequency-to-voltage  converters  or  other 
devices  which  require  a  programmable- frequency  square  wave 
input . 

Figure  2  [Ref.  2]  shows  a  detailed  block  diagram  of 
the  timer  pacer  card.   In  the  Light  Scattering  experiment 
the  timer  pacer  card  is  programmed  to  the  continuous  mode 
and  drives  the  external  trigger  input  of  2  high  speed  A/D 
converter  cards.   The  timer  pacer  is  itself  externally  trig- 
gered by  the  Reticon  Photodiode  array  #1  which  also  syn- 
chronizes the  operation  of  Photodiode  array  #2. 
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2 .   Analog  to  Digital  Converter  Card  #69751A 

The  69751A  is  a  12  bit,  analog  to  digital  (A/D)  con- 
verter used  to  measure  bipolar  DC  voltages  in  one  of  four 
ranges:   ±100  MV,  ±1  V,  ±10  V,  or  ±100  V.   A  voltage  measure- 
ment can  be  initiated  by  a  programmed  instruction  or  by  an 
external  trigger  signal.   Three  manual  range  switches  on  the 
card  select  either  the  ±100  MV,  ±1  V,  or  the  ±10  V  range. 
Input  voltages  in  the  ±100  V  range  are  connected  to  the  card's 
divide-by-ten  attenuator.   Sixteen  optional  DC  input  ranges 
(8  bipolar  and  8  unipolar)  are  also  available.   Instructions 
for  modifying  and  recalibrating  the  A/D  card  for  the  optional 
DC  input  ranges  are  described  in  section  11  of  Ref.  2. 

The  A/D  card  uses  a  high  performance  sample  and  hold 
amplifier  and  a  successive  approximation  A/D  converter  which 
provide  a  high  conversion  speed  and  excellent  accuracy  charac- 
teristics.  The  DC  input  voltage  is  guarded  to  maintain  high 
input  impedance.   Optically  coupled  isolators  are  used  to 
isolate  the  DC  input  voltage  from  the  data  lines. 

The  12  bit  digital  word,  indicating  the  magnitude  and 
sign  of  the  measured  voltage,  is  stored  on  the  card  and  is 
always  available  for  readback  to  the  multiprogrammer  main- 
frame memory  using  card  subaddress  0.   As  many  as  1339  read- 
ings can  be  stored  at  one  time  in  mainframe  memory.   The 
multiprogrammer  firmware  converts  the  digital  word  decimal 
form  for  readback  to  the  controller.   The  controller  can 
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read  back  the  voltage  value  using  appropriate  HP-IB  extended 
talk  address. 

The  use  of  Extended  Talk  Addresses  and  other  second- 
ary addressing  with  the  HP-85  computer  requires  the  custom 
HP-IB  command  "SEND"  which  is  described  in  detail  in  the 
HP-IB  I/O  programming  manual  [Ref.  3]. 

The  12  bit  data  word  (voltage  reading)  is  routed  to 
the  A/D's  edge  connector  so  that  the  A/D  can  be  used  in  con- 
junction with  the  69790B  memory  card  to  make  measurements  at 
rates  that  exceed  the  throughput  of  the  controller.   Up  to 
33000  readings  per  second  can  be  transferred  from  the  A/D 
converter  to  the  memory  card  via  the  edge  connectors  on  each 
card.   Each  reading  can  be  initiated  by  program  commands  or 
by  an  External  Trigger  (EXT)  signal. 

Either  the  69736A  Timer  Pacer  or  the  69735A  Pulse 
Train  Output  card  can  be  used  to  externally  trigger  the  A/D 
converter.   The  69736A  Timer/Pacer  is  programmed  in  a  con- 
venient seconds-milliseconds-microseconds  format  for  pro- 
grammable delay  or  continuous  triggering. 

Figure  3  shows  the  functional  block  diagram  of  the 
69751A  Analog  to  Digital  converter  card.   Various  subad- 
dresses  can  also  be  written  using  the  software  described  in 
section  5  of  the  User's  Guide  [Ref.  2]. 
3.   Memory  Card  #6  9790B 

The  69790B  Memory  Card  is  a  FIFO  (first  in  first  out) 
buffer  that  can  handle  a  continuous  stream  of  data  that  is 
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being  transferred  from  the  outside  world  to  the  controller 
or  from  the  controller  to  the  outside  world.   In  another 
operating  mode,  the  recirculating  input  mode,  it  can  act  as 
a  circular  buffer  that  continuously  stores  data  coming  in 
from  the  outside  world  until  it  is  programmed  to  stop  ac- 
cepting new  data.   At  this  time  it  contains  the  last  'n' 
words  of  data  received,  'n'  being  the  capacity  of  the  card. 
One  more  operating  mode,  the  recirculating  output  mode,  al- 
lows all  or  part  of  the  memory  contents  previously  acquired 
to  be  read  repeatedly  by  an  external  device. 

Each  69790B  memory  card  consists  of  two  plug-in 
cards  that  occupy  two  adjacent  slots  in  the  mainframe.   The 
card  designated  card  #1  occupies  the  left-hand  slot  and  com- 
municates with  card  #2  through  a  cable. 

A  memory  card  assembly  can  communicate  bidirection- 
ally  with  the  controller  and  with  an  external  device  con- 
nected to  the  edge  connector  on  card  #1  of  the  pair.   External 
data  transfers  through  16  input  or  16  output  lines  and  are 
controlled  by  means  of  three  input  or  three  output  handshake 
lines.   Special  purpose  Memory  Input  (MI)  and  Memory  Output 
(MO)  instructions  allow  data  transfers  to  be  made  between 
the  controller  and  a  memory  card  assembly  without  storing 
the  data  in  mainframe  memory. 

In  a  system  that  contains  more  than  one  memory  card 
assembly,  each  assembly  can  handle  data  transfers  to  or  from 
the  outside  world  simultaneously  with  external  data  transfers 
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being  made  by  the  other  assemblies.   Absolutely  no  inter- 
actions occur  between  their  external  interfaces.   An  external 
interface  is  independent  not  only  of  the  other  external  in- 
terfaces, but  also  of  any  other  card  or  mainframe  activity. 
Thus,  within  a  single  assembly,  an  external  device  can  write 
to  a  memory  card  while  the  controller  is  reading  from  another 
part  of  the  same  card,  or  an  external  device  can  read  from  a 
memory  card  while  the  controller  is  writing  to  another  part 
of  the  same  card.   Through  the  use  of  controller  interrupts 
or  another  means  of  program  control,  data  can  be  loaded  into 
memory  and  read  from  memory  on  a  continuous  basis  in  either 
an  input  or  an  output  operation. 

Figure  4  is  a  functional  block  diagram  of  the  69790B 
memory  card. 
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Figure  4.   Multiprogrammer — Memory  Board  [Ref.  2] 
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III.   PARTICLE  SIZE  DETERMINATION  BY  LIGHT  SCATTERING 
OF  A  HELIUM  NEON  LASER  BEAM 

The  purpose  of  this  experiment  is  to  determine  the  volume 
to  surface  mean  diameters  (D-32)  ^f  particulate  matter  in  the 
exhaust  plume  and  combustion  chamber  of  a  solid  propellant 
rocket  motor. 

The  equipment  test  set  up  is  shown  in  Figure  5.   Refer- 
ence 4  provides  extensive  information  on  the  physical  princi- 
ples involved  and  a  description  of  the  optics  system. 

The  helium  neon  laser  is  remotely  controlled  under  soft- 
ware control  and  radiates  a  beam  at  a  wavelength  of  0.6  32  8 
micrometers.   The  beam  is  then  expanded  and  columated.   It 
is  then  passed  through  a  limiting  orifice  to  reduce  the  beam 
diameter  to  2.2  millimeters. 

The  beam  is  subsequently  passed  into  a  beam  splitter 
where  one  beam  is  turned  through  90°  and  routed  through  a 
window  in  the  rocket  motor  casing.   The  other  beam  is  aligned 
directly  in  the  exhaust  stream  of  the  rocket  motor  approxi- 
mately two  inches  behind  the  nozzle  throat. 

Each  beam  is  passed  through  a  focusing  lens  and  directed 
through  a  narrow  pass  filter  to  the  edge  of  a  linear  array 
of  photodiodes. 

Two  EG&G  RETICON  G  Series  Solid  State  Line  Scanners  are 
used,  each  containing  1024  photodiode  elements.   The  Reticon 
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RL1024G  series  photodiode  array  is  optimized  for  solid  state 
image  sensor  applications.   It  contains  1024  diodes  arranged 
in  a  linear  array  on  25  micron  centers.   Diode  #1  is  normally 
set  1.54  millimeters  off  of  the  laser  beam  centerline  to  re- 
duce the  amount  of  transmitted  light  present,  and  to  allow 
for  increased  sensitivity  to  the  diffracted  light  pattern. 

The  G  series  array  is  connected  through  an  RC-100/106 
sample  and  hold  "Boxcar"  circuit  which  provides  its  own 
clock  oscillator  for  self  scanning,  although  provisions  for 
externally  pacing  the  array  with  an  external  clock  signal  are 
installed. 

Photodiode  array  #1  (Fig.  5)  is  operated  in  the  internal 
clock/internal  start  mode  and  is  completely  independent.   The 
internal  clock  oscillator  is  set  to  approximately  30  KHz.   To 
prevent  integrated  "dark  current"  from  making  a  significant 
contribution  to  the  output  charge  the  total  time  between 
start  pulses  should  be  kept  less  than  about  40  milliseconds. 

The  countdown  circuit  on  the  array  board  counts  'n' 
clock  pulses  to  generate  a  start  pulse  in  the  internal  start 
mode.   The  countdown  switches  are  initially  set  according  to 
the  formula: 

n  =  t^  •  f 
L     c 

where  t   =  time  between  start  pulses  in  seconds 

Li 

f   =  clock  frequency  in  Hz. 
n  =  switch  setting 
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e.g.    n  =  36  (lO"^)  •  30  (10^) 
n  =  1080 
A  countdown  switch  setting  of  1080  to  1  will  provide  a  time 
between  start  pulses  of  36  milliseconds  at  a  clock  frequency 
of  30  KHz,  and  excessive  "dark  current"  error  build-up  is 
avoided. 

Photodiode  array  #2  monitors  the  second  laser  beam  which 
emerges  from  a  window  in  the  rocket  motor  body.   Array  #2 
operates  in  the  external  clock/external  start  mode  and  is  a 
slave  to  photodiode  array  #1.   Therefore,  the  countdown  cir- 
cuits in  array  #2  are  not  used.   Array  #2  has  only  one  out- 
put, the  video  "boxcar"  signal  representing  the  photodiode 
voltage  output  which  is  routed  directly  to  A/D  converter  #2 
in  the  6942A  multiprogrammer . 

The  boxcar  voltage  outputs  from  the  photodiode  arrays 
are  routed  through  analog  low  pass  filters  which  pass  all 
signal  components  less  than  3000  Hz.   20  db  gain  is  attained 
in  the  low  pass  filter  amplifier.   The  filtered  outputs  are 
routed  to  the  voltage  input  terminals  {+W,  -Y)  of  the  high 
speed  analog  to  digital  converters  located  in  the  6942A 
multiprogrammer.   When  a  series  of  readings  are  desired  the 
HP-85  programs  the  69736A  timer  to  output  a  square  wave  output 
at  approximately  30  KHz  when  externally  triggered  by  the 
blanking  pulse  output  from  photodiode  array  #1. 

The  square  wave  output  from  the  timer  pacer  card  is  di- 
rected to  the  external  trigger  inputs  of  both  69751A  high 
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speed  (A/D)  cards  which  begin  making  1024  conversions  of  the 
analog  voltage  present  at  terminals  (+W,  and  -Y) .   The  12  bit 
digital  outputs  from  terminals  A  through  N  are  hard  wired  to 
the  corresponding  edge  connector  pins  of  the  69790B,  4K  memo- 
ry  cards  associated  with  each  (A/D)  card.   The  memory  is  pro- 
grammed to  accept  a  12  bit  format  and  to  store  1024  readings 
in  a  first  in  first  out  (FIFO)  input  mode.   The  memory  card 
is  then  armed  to  interrupt  the  computer  after  1024  readings 
have  been  stored  on  the  card.   The  (A/D)  card  performs  a  con- 
version of  the  analog  voltage  at  its  input  terminals  (W,  Y) 
every  time  it  observes  a  trigger  from  the  timer  pacer  card  at 
its  external  trigger  terminals  (m)  at  a  rate  of  30  KHz.   The 
(A/D)  card  signals  the  completion  of  each  reading  by  output- 


ting  an  end  of  conversion  (EOC)  and  its  complement  (EOC) . 


The  EOC  is  connected  to  the  data  available  (DAV)  terminal  of 
each  69790B  memory  board. 

The  procedure  employed  for  an  experiment  is  outlined  be- 
low.  Figure  6  presents  a  flow  diagram  of  the  system.   After 
applying  power,  the  HP-85  is  loaded  with  a  program  disc  in 
drive  #0  and  a  data  disc  in  drive  #1.   The  program  named 
"HELP"  is  loaded  and  run  to  initialize  variables  and  provide 
documentation.   At  the  completion  of  "HELP"  the  next  program 
is  automatically  loaded  and  executed  by  a  process  called 
chaining,  which  allows  a  program  of  virtually  unlimited  length 
to  be  run  in  the  limited  memory  space  of  32K  of  the  HP-8  5 
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Figure    6.      Light    Scattering   Program  Flow   Chart — Part   One 
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computer.   "ACQDTA"  now  performs  the  process  control,  system 
timing,  and  data  acquisition  tasks. 

First,  the  laser  is  energized,  data  files  are  created 
with  user  specified  file  names,  and  initial  readings  are 
taken  from  both  photodiode  arrays  by  software  command.   Each 
array  simultaneously  sends  its  output  to  its  69751A  (A/D) 
converter  in  the  6942A  multiprogrammer .   Each  (A/D)  card  is 
externally  triggered  every  33  microseconds  by  the  69736A 
Timer/Pacer  card  and  converts  the  analog  voltage  present 
at  terminals  (+W,  -Y) .   The  resulting  12  bit  digital  out- 
put is  sent  to  a  69790B  (4K)  memory  card.   Each  memory  card 
is  programmed  to  store  sequentially  1024  readings  and  to 
interrupt  the  computer  after  1024  readings  have  been  stored. 
At  this  time  the  numerical  data  are  transferred  via  the 
HP-IB  to  the  controller,  first  array  #1  then  array  #2.   This 
raw  data  represents  the  "no  particle"  or  "zero  scattering" 
condition. 

The  "no  particle"  numerical  arrays  are  then  stored  as 
a  data  file  on  the  disc  to  preclude  loss  of  the  data. 

The  computer  then  performs  a  continuity  check  of  the 
firing  circuit  using  the  3456A  in  a  resistance  measurement 
mode.   After  a  successful  continuity  check,  the  3456A  is  pro- 
gramm.ed  to  record  DC  voltage.   It  is  connected  to  the  firing 
switch  while  displaying  (on  the  HP-8  5) ,  "STANDING  BY  FOR 
IGNITION" ,  and  energizes  the  nitrogen  purge  for  the  motor 
windows.   When  the  operator  closes  the  ignition  switch  a 
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voltage  greater  than  10  volts  is  sensed  by  the  3455A  and 
causes  the  computer  to  energize  the  drive  motor  of  the  visi- 
corder,  which  records  the  chamber  pressure-time  trace. 

Once  the  computer  senses  the  ignition  sequence  it  con- 
nects the  3456A  to  the  pressure  transducer  in  parallel  with 
the  visicorder  trace. 

When  the  combustion  pressure  rises  to  a  user  specified 
threshold  pressure  in  (PSI)  the  computer  delays  or  waits  a 
user  specified  number  of  seconds,  and  triggers  the  multipro- 
grammer  to  acquire  a  new  data  set  from  the  photodiode  arrays. 
When  the  data  are  taken  a  pulse  output  voltage  (TTL)  is  sent 
to  another  visicorder  channel  as  an  event  marker  signal  to 
indicate  data  acquisition  time  on  the  pressure-time  trace. 
Readings  taken  during  the  rocket  motor  burn  time  will  be 
different  than  pretest  readings  since  particulate  matter  is 
present  in  the  gas  stream.   This  data  is  transferred  to  the 
computer  and  then  to  a  data  file  on  the  data  disc  as  before. 

After  5  seconds  the  visicorder  is  disengaged,  the  window 
nitrogen  purge  turned  off,  and  the  laser  is  deenergized.   The 
next  program  is  chained. 

"CMBINE"  is  automatically  loaded  and  run.   This  program 
executes  disc  operations  which  read  the  "no  particle"  data 
and  the  "particle"  data  into  the  computer.   It  then  separates 
this  data  into  photodiode  array  #1  data  and  photodiode  array 
#2  data  with  the  same  file  names  as  previously  specified  by 
the  user.   Each  file  contains  1024  records,  each  containing 
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two  numbers,  a  "no  particle"  and  a  "particle"  voltage.   The 
computer  prints  a  reminder  of  which  data  is  in  which  data 
file  and  then  chains  the  next  program  segment. 

"PLTDTA"  plots  the  data  from  a  data  file  specified  by 
the  user,  plotting  diode  voltage  output  versus  diode  number 
(Fig.  7)  .   The  operator  can  elect  to  plot  the  other  data 
file  or  load  and  run  the  data  reduction  program. 

The  operator  is  then  prompted  to  load  and  run  "REDUCE"  to 
begin  the  reduction  of  the  raw  data.   Figure  8  shows  a  flow 
chart  for  the  data  reduction  process.   "Reduce"  functions  to 
compute  a  normalized  intensity  vector  and  an  angle  theta 
vector.   The  theta  vector  is  computed  based  on  diode  number 
(1  to  1024)  ,  lens  focal  length,  and  the  index  of  refraction 
of  the  media  containing  the  particles  by  the  formula: 

K^  +  n( .025) 

9  =  = radians 

F  m 

where  K^  =  the  distance  from  the  centerline  of  the  laser  beam 
to  the  centerline  of  diode  #1  (mm) 

n  =  diode  number  (1  to  1024) 

F  =  lens  focal  length  (mm) 

m  =  index  of  refraction 
air  =1.0 
water  =1.35 

.025  =  distance  between  diode  centers  (mm) 
Upon  examination  of  the  "PLTDTA"  results,  the  location 
of  the  peak  intensities  are  printed.   The  user  is  then  prompt- 
ed for  the  diode  locations  where  the  slope  is  to  be  determined 
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and  then  projected  back  to  obtain  a  centerline  value  for  the 
scattered  laser  power  using  a  linear  regression  technique. 

First,  the  "zero  scattering"  data  is  subtracted  from  the 
"with  particle"  data.   The  resulting  1024  element  vector  is 
then  regressed  against  diode  number  to  obtain  the  intensity- 
axis  intercept  representing  the  centerline  value  by  the  fol- 
lowing formula: 

nA  +  BX^  =  Y^ 

X^A  +  BX2  =  Y2 

where   A  =  intensity-axis  intercept 

n  =  number  of  data  points 

B  =  slope  of  intensity  vs  theta  line 

n 

X,  =   y  0. 
1    •  ^-,  1 


i=l 
n 
X 

i=l 
n 
Y^ 

i=l 
n 
Y^ 

i=l 


2    ■  ^-,  1 


1  =  .^.^i 


2    .^,11 


and  subsequently,  a  normalized  intensity  vector  is  formed  as 

■^9N0RiM      ,„       „     . 
(^CLp  -  ^CLnp) 

where  i=ltol024 

Ifl^-r^T,..  =  normalized  intensities 
"NORM 

Vg  (i)  =  voltage  with  particles  (vector) 
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Vqj^  (i)  =  voltage  without  particles  (vector) 
(VcLp  -  '^CLnp)  "  ^  (previous  calculation) 

The  normalized  intensity  is  then  plotted  on  the  7225B 
plotter  versus  the  theta  vector  in  radians  (Fig.  9) .   This 
data  is  then  stored  on  disc  for  future  use  in  a  user  speci- 
fied file.   The  operator  is  then  asked  if  this  run  is  to  be 
a  calibration  curve;  if  yes,  then  a  curve  fitting  program 
"CALIB"  is  automatically  loaded  and  executed. 

"CALIB"  performs  a  tenth  order  polynomial  curve  fit  to 
the  data  set  specified  by  the  user  in  order  to  generate  a 
smooth  curve  for  future  comparison  purposes  with  both  the 
universal  (i.e.  theoretical  [Ref.4])  curve  ("UNICUR")  and 
other  raw  data  runs.   The  user  first  enters  the  "D  bar" 
value,  i.e.  the  average  diameter  of  the  particles  being  used 
for  the  calibration. 

A  theta  bar  vector  is  then  computed  from  the  theta  vector 
as  follows:  , 

-'       I' 

i 

A 

where   6.  =  theta  vector  (1-1024)  (radians) 

D-.-P  =  D  bar  (microns)  mean  volume  to  surface  diameter 
A  =  .6328  (microns)  wavelength 
Normalized  intensity  is  then  plotted  versus  the  non- 
dimensional  theta  bar  (6)  as  in  figure  10.   The  user  is  then 
asked  to  enter  the  minimum  and  the  maximum  theta  bar  of 
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interest.   The  data  is  regressed  to  form  a  smooth  curve  fit 
with  a  polynomial  of  up  to  tenth  order  as  specified  by  the 
user  and  taking  the  form: 

I^,  =  b-  +  b^  9  +  b^9^  +  b^9^  +  .  ,  ,  +  b  o"^ 
N     0     1      2       3  m 

where  the  coefficients  are  obtained  by  solving  the  system  of 
m+1  equations 

nb^  +  b^  y  9,  +  b^  y  9.   +  .  .  .  +  b   y  9.^^  =   I  !• 
^     ■^i=l  ^     ^i=l  ^  "^i=l  "■  i=l  ^ 

n         n^nS  ^^mj-i^i 

b.  I  9.  +  b,  I  9/  +  b^  I  9.   +  .  .  .  +  b^  y  9.""^^  =   y  1.9. 


1-  V  n  ni  ,  I-  "  ^m+1  ,  ,  ,  VA2m   r^.  o,m 

b-)9.  +   b^)^  +  .     .     .    +b)9    =)I.9 

O^il"  m^       ^1 

Once  the  polynomial  coefficients  have  been  determined,  a 
smooth  100  element  curve  is  computed,  plotted  and  stored  on 
the  disc  in  a  user  specified  file. 

The  smooth  curve  is  plotted  on  top  of  the  intensity  vs 
theta  bar  plot  (Fig.  10)  to  show  the  correlation,  which  is 
routinely  excellent.   If  desired,  this  data  file  containing 
smooth  curve  data  can  be  saved  permanently  on  the  floppy  disc 

"REDUC2"  is  automatically  loaded  and  executed  upon  com- 
pletion of  "CALIB".   The  function  of  "REDUC2"  is  to  plot  the 
universal  theoretical  curve  (datafile  "UNICUR")  on  a  three 
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cycle  semi-logarithmic  scale  (Fig,  11) .  After  the  theoretical 
curve  is  plotted  the  user  may  elect  to  plot  any  other  calibra- 
tion curve  for  comparison  purposes  (Fig.  11). 

At  the  conclusion  of  "REDUCE",  if  the  operator  had  not 
designated  a  run  as  a  calibration  run,  the  next  data  reduc- 
tion program  is  chained. 

"REDUC3"  functions  to  compute  D^-  f^om  the  normalized 
intensity  and  theta  data  by  comparing  the  data  to  specified 
reference  curves  (UNICUR,  etc.)  at  15  preset  intensity  levels. 
An  average  D^-  is  computed  across  this  range  of  intensity 
levels.   REDUC3  has  tabular  output  only. 

REDUC4  can  be  run,  if  desired,  to  plot  a  data  run  on  a 
semi-log  scale  with  a  user  defined  D^»  to  compare  with  the 
Universal  curve. 
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IV.   VITIATED  AIR  HEATER 

A  schematic  of  the  vitiated  air  heater  system  at  the  com- 
bustion laboratory  is  shown  in  figure  #12.   It  was  designed 
to  supply  air  at  any  desired  mass  flow  rate  within  the  system 
limitations  (approximately  4  Ibm/sec)  and  at  any  temperature 
up  to  approximately  1500 °R.   The  heater  output  air  is  routed 
to  the  input  sections  of  the  Solid  Fuel  Ramjet  or  the  Turbo- 
jet combustor. 

The  vitiated  air  heater  is  fueled  by  gaseous  ethylene  and 
bottled  oxygen  replaces  that  consumed  in  the  combustion  pro- 
cess.  Large  high  pressure  air  tanks  act  as  the  primary  air 
supply  and  the  mass  flow  rate  is  determined  using  a  sonically 
choked  converging  nozzle. 

The  data  acquisition  system  was  designed  to  monitor  and 
control  the  operation  of  the  air  heater  in  order  to  provide 
the  desired  air  flow  rate  and  temperature  and  to  provide 
built-in  system  safety. 

Two  Hastings  precision  flowmeters  and  two  Hastings  flow 
controllers  were  installed  to  accurately  meter  and  control 
the  oxygen  and  ethylene  flow  rates  to  achieve  the  desired 
fuel-air  ratio  and  oxygen  make-up  flow  rate.   Figure  #13 
shows  the  interconnection  of  the  flowmeters  and  controllers. 

A  description  of  the  system  operation  is  presented  below. 
The  'HEATER'  program  disc  is  loaded  into  the  HP-85  data 
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Figure    13.      Flow  Control   Block   Diagram 
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acquisition  system.   The  "SET  UP"  option  is  selected  to  ini- 
tialize variables  and  perform  preliminary  calculations  to 
determine  initial  manual  control  settings.   Figure  14  depicts 
the  flow  chart  for  the  preliminary  portion  of  the  'HEATER' 
program. 

The  operator  inputs  the  desired  mass  flow  rate  of  air  in 
lbs  mass/sec.   Other  inputs  are  high  pressure  air  tank  pres- 
sure (0-3500  PSIG) ,  the  estimated  cold  air  temperature  (°F) , 
desired  hot  air  temperature  (°R),  and  ambient  pressure  (inches 
Hg)  . 

The  required  fuel/air  ratio  is  then  computed  for  ethylene 
and  air  from  a  regression  of  data  obtained  from  the  "PEPCODE" 
chemical  equilibrium/adiabatic  combustion  program  [Ref.  5]. 

p    ^HOT  -  "°    *FUEL 
78000    -  ri,^^j, 

Once  F  is  obtained,  the  ethylene  and  oxygen  flow  rates 
are  given  by: 

m^  u   =  F  •  m_    (#  mass /sec) 
^2^4        ^2 

For  the  complete  combustion  of  ethylene  with  oxygen  in  air: 

C2H^  +  3O2  ^  2CO2  +  2H2O 

Thus,  the  desired  mass  flow  ratio  of  oxygen  to  ethylene  is 

m   =  3.429  •  m 
O2  C2H^ 
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The  volume  flow  rates  of  oxygen  and  ethylene  in  SCFM  are 

Q    =  734.89  •  m^ 

Q      =  839.87  •  m 
^2^4  ^2^4 

thus  Qp.   =  3Q 

^2      ^2^4 

where  Q  is  the  volume  flow  rate.   The  oxygen  and  ethylene 
flowmeters  produce  0-5  Volt  DC  signals  linear  with  volume 
flow  rate  (0-20  SCFM  for  oxygen  and  0-6.5  SCFM  for  ethylene) 
The  required  oxygen  and  ethylene  set  voltages  are  computed 
from 

V^   =  183723  •  m    (Ibm/sec)  (millivolts) 
^2  ^2 

V^  „   =  646054  •  m^  „   (Ibm/sec)  (millivolts) 
^2^4  ^2^4 


V   is  generated  by  the  dual  channel  digital  to  analog 
^2 
(D/A)  converter  card  in  the  3497A  and  sent  to  the  REMOTE 

terminals  of  the  oxygen  flow  controller  (Fig.  13) .   This  flow 

controller  is  set  to  the  'BIAS'  mode  of  operation  with  the 

Local  Set  Point  (LSP)  control  set  to  zero. 

In  this  position  (Fig.  15a)  the  input  stages  of  the  flow 

controller  sum  the  three  voltages  and  generate  an  error 

signal: 

ERROR  =  INPUT  -  LSP-  REMOTE 

This  signal  drives  the  flow  control  until  the  error  signal  is 
equal  to  zero,  which  occurs  when  the  desired  oxygen  flow  rate 
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is  achieved.   With  the  LSP  set  to  000,  the  oxygen  set  voltage 
from  the  D/A  converter  acts  as  a  remote  set  point  to  control 
the  desired  flow  rate.   The  upper  section  of  figure  13  shows 
the  interconnections  of  the  second  flow  controller,  the  ethyl- 
ene flowmeter,  and  the  oxygen  flowmeter  in  the  RATIO  mode  of 
operation.   In  this  mode,  the  flow  controller  acts  to  provide 
the  proper  ratio  of  ethylene  to  oxygen  under  all  conditions 
by  controlling  the  ethylene  volume  flow  rate  to  maintain  a 
1:3  ratio.   Thus,  the  computer  controls  the  oxygen  flow  rate 
which  in  turn  controls  the  ethylene  flow  rate. 

Figure  15b  shows  the  flow  controller  input  circuit  in  the 
1-5  RATIO  mode.   In  this  case,  the  REMOTE  (OXYGEN)  signal  is 
equal  to  the  product  of  the  Local  Set  Point  (LSP)  and  the 
INPUT  (ETHYLENE)  signal. 

LSP  •  INPUT  =  REMOTE 

However,  the  input  (ethylene)  and  remote  (oxygen)  flow 
capacities  are  different.  Dividing  by  the  full  scale  flow 
capacities  in  SCFM  normalizes  the  equation. 

^  ^^    /  INPUT  \    / REMOTE  \ 

REMOTE  =  Oxygen  volume  flow  rate 
INPUT  =  Ethylene  volume  flow  rate 


Then, 

INPUT 
REMOTE 

Thus,  LSP  =  0.975. 


\  20  /VLSP/    3 
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For  example,  if  the  desired  oxygen  flow  rate  is  3.5  SCFM 
or  0.00476  Ibm/sec,  then  the  OXYGEN  set  voltage  is 

V^   =  3.5  SCFM  L^c^^M   =  0.875  Volts 
0^  20  SCFM 

The  computer  sends  0.875  volts  to  the  remote  terminals  of 
the  oxygen  flow  controller.   The  RATIO  controller  responds  by 
driving  the  flow  control  valve  until  the  ethylene  flowmeter 
reads: 

INPUT  •  LSP  =  REMOTE 

INPUT  =  ?,'lll    =  0.8974  VOLTS 

This  corresponds  to  an  ethylene  volume  flow  rate  of 

Q^  ^   =  0.8974  Volts  ^.'  I,  ^^J^ 

C2H.  3  Volts 

Q_  „   =1.166  SCFM 
^2^4 

which  is  1/3  of  the  oxygen  volume  flow  rate. 

This  completes  the  preliminary  calculations  and  the  sys- 
tem is  now  turned  on  to  provide  air  flow  without  any  heat 
addition.   Thus,  the  operator  selects  the  'FLOW  option  (Fig. 
14)  and  the  computer  measures  ambient  pressure  then  opens  the 
air  valve,  allowing  air  to  flow  to  the  heater  while  other 
components  are  maintained  in  the  off  position.   High  pressure 
air  passes  through  a  choked  flow  nozzle,  where  P-p-,  and  T-p-i 
are  measured.   Based  on  these  measurements  the  actual  mass 
flow  rate  is  computed  from: 
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mj^-j-j^  =  ,  —         (Ibm/sec) 


R   T 


Tl 


^c 


where 


r 

Y+1 

^1 

= 

kt.r 

Y 

= 

1.4 

AIR 

A* 

= 

TTCl 

4 

(in^) 

d* 

= 

.235 

(in)    typical 

^D 

= 

0.97 

,    discharge   coef. 

R 

53-3 

ft-lbF 

■"'•■'      lbm°R 
Tr£i      =  measured   temperature,    °R 

Prp-  =  measured  pressure,  PSIA 

--^0    0      IbF  ft 

^c   -  32.2  2 

Ibm  sec 

The  actual  mass  flow  rate  of  air  is  displayed  and  an  up- 
dated oxygen  set  voltage  is  computed  and  set.   The  required 
pressure  correction  AP<p.  is  computed  to  meet  the  desired 
mass  flow  rate. 

The  operator  may  then  accept  the  flow  conditions  and 
proceed  to  another  program  segment.   If  he  accepts  these 
flow  conditions,  the  actual  mass  flow  rate  of  air  is  set 
equal  to  the  desired  mass  flow  rate.   If  he  is  not  satisfied 
with  existing  flow  conditions,  P-p   is  adjusted  manually  by 
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applying  the  recommended  AP-p  .   The  'FLOW  option  is  auto- 
matically executed  again  until  the  desired  accuracy  is 
attained. 

Once  the  desired  mass  flow  rate  of  air  is  obtained,  the 
'START'  option  is  executed.   Figure  16  shows  a  flow  chart  of 
this  program  segment.   This  routine  first  closes  all  valves 
and  measures  P>p  which  represents  ambient  pressure  under  no 
flow  conditions.   The  proper  oxygen  set  voltage  is  set  ac- 
cording to  the  latest  flow  calculations.   The  computer  then 
turns  on  the  high  pressure  air  flow,  energizes  the  igniter 
system,  and  turns  on  the  ethylene  and  oxygen  flow  to  the 
heater  section.   The  output  hot  air  temperature  is  monitored 
by  a  thermocouple  and  measured  by  the  computer  at  0.5  second 
intervals  during  the  start  up  sequence,  beginning  3  seconds 
after  ignition.   The  igniter  is  turned  off  after  5  seconds  of 
operation  unless  secured  earlier  by  manual  control. 

During  heater  start  up  a  typical  temperature-time  trace 
would  be  as  shown  in  figure  17. 

Temperature  is  closely  monitored  during  the  start  se- 
quence and  igniter  shutdown  is  executed: 

1)  when  Tj^Qrp  j^j-^   reaches  700°R  (if  not  previously  secured 
by  the  5  second  timer) 

2)  if  Thot  air  fslls  by  50°R  in  any  5  second  interval  then 
HEATER  shutdown  is  called  to  close  the  OXYGEN  and  ETHYL- 
ENE valves,  close  igniter  valves,  set  O2  set  voltage  =  0 
and  close  the  high  pressure  air  valve. 

^^^^  ^HOT  AIR  ^^^^hes  T^^^^^^^    -  100°  the  Heater  Mainten- 
ance  routine  is  called,  terminating  the  start  sequence.   The 
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Heater  Maintenance  Program  is  assigned  to  timer  #3  and  set  to 
automatically  execute  every  20  seconds,  while  another  project 
program  segment  is  run. 

Figure  18  shows  the  flow  chart  for  the  Heater  Maintenance 
Program  segment  (HTRMNT) ,  which  is  assigned  to  Timer  #3  to 
interrupt  the  computer  at  20  second  intervals  to  control  the 
heater  output  temperature.   On  interruption,  the  computer 
measures  the  output  air  temperature  and  compares  it  to  the 
desired  air  temperature.   If  the  difference  is  greater  than 
200 °R  the  heater  shutdown  program  is  called  and  secures  the 
heater  entirely.   If  the  difference  is  less  than  100°R  the 
computer  records  the  air  temperature  and  returns  to  whatever 
routine  was  operating  prior  to  the  interrupt.   However,  if 
the  difference  between  actual  temperature  and  desired  tempera- 
ture is  over  100°  and  less  than  200°R  the  sign  of  the  error 
is  determined  and  the  fuel  air  ratio  is  adjusted  to  compen- 
sate for  the  temperature  deviation.   Oxygen  and  ethylene  flow 
rates  are  updated  and  a  new  oxygen  set  voltage  is  sent  to  the 
oxygen  flow  controller  to  effect  a  change  in  the  output  tem- 
perature.  The  output  temperature  is  then  recorded  and  a  re- 
turn from  interrupt  is  executed. 
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INTERRUPT 
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HTRMNT 


1)  Measure  T^^^  ^^^ 

2)  Compute  AT=T^^-TJ^^^^ 


YES 


YES 


HEATER 
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TOO  COOL 


5 


TOO  HOT 


T)  I^F=     T/78000 

2)  F  =  F+  AF 

3)  Update  otq^  +  ^2^4 

4)  Vo2  =  13373  •  iho^ 

5)  Update  0^  set  voltage 


3 


Record  T 


HOT  AIR 
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Figure  18.   HEATER  Flow  Chart — Part  Three 


61 


V.   SOLID  FUEL  RAMJET 

The  Solid  Fuel  Ramjet  (SFRJ)  test  facility  at  the  combus- 
tion laboratory  has  been  constructed  to  facilitate  the  mea- 
surement of  performance  as  affected  by  design  and  operating 
variables.   Figure  19  presents  a  schematic  of  the  SFRJ  motor 
which  is  mounted  on  a  thrust  stand.   Figure  20  presents  a 
schematic  of  the  air  supply  system.   Air  flow  is  directed  to 
the  head-end  and  sometimes  also  to  the  aft  mixing  chamber 
(bypass-air)  of  the  combustor.   In  general,  a  test  sequence 
is  as  follows:   air  flow  rate(s)  are  turned  on,  the  igniter 
is  actuated,  a  small  amount  of  ethylene  is  mixed  with  the 
head-end  air  flow,  the  motor  ignites,  the  igniter  and  ethyl- 
ene flows  are  terminated,  the  motor  runs  for  10-45  seconds 
duration,  air  flow  is  terminated  and  nitrogen  gas  is  flowed 
through  the  fuel  grain  for  five  seconds  to  ensure  extinguish- 
ment of  the  flame.   During  the  test  it  is  important  to  mea- 
sure air  flow  rates,  inlet  air  temperature,  head-end,  and 
mixing  chamber  pressures.   Postfire  data  collected  include 
fuel  grain  weight  loss  which  is  needed  to  calculate  the 
average  fuel  flow  rate  and  average  fuel  regression  rate. 

Temperature  rise  combustion  efficiency  for  the  test  is 
determined  from 

T+-     -  Tt  . 
_   ^exp "-air 

'^AT    Tt^u  -  Tt  . 
^th     ^air 
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where  T^.     is  the  measured  inlet  air  stagnation  temperature, 
Ttj-h    ^^   ^^^  theoretical  combustion  temperature  determined 
from  adiabatic,  equilibrium  combustion  calculations  using 
the  measured  fuel-air  ratio,  inlet  air  temperature  and  aver- 
age chamber  pressure.   Tt     is  the  experimentally  determined 
combustion  temperature.   The  latter  is  calculated  from  con- 
servation of  mass,  using  gas  properties  determined  from  the 
equilibrium  calculations.   Thus, 


-exp 


S^c\h 


-|2 


m, 


R 


ilk) 


Y-l 


where   m^  =  m^.^  +  m^^^^ 

A^,  =  nozzle  throat  area 
th 

P_  =  average  chamber  stagnation  pressure 

C-.  =  nozzle  discharge  coefficient 

Average  fuel  regression  rate  is  calculated  based  upon 
fuel  weight  loss  during  the  burn.   Thus, 


-4 


4AW  ^  a  2 

irLp     1 


and 


r  = 


f 1 

2t^ 


where  AW  =  fuel  weight  loss 
L  =  fuel  grain  length 
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p  =  fuel  density 

d^  =  final  average  fuel  internal  diameter 

J. 

t,  =  burn  time 

D 

f  =  average  fuel  regression  rate 

The  average  chamber  pressure  (P  )  used  to  determine  T-j- 
is  obtained  from  the  pressure  time  trace  recorded  during  the 
test.   A  typical  trace  is  shown  in  figure  21. 

The  test  control/data  acquisition/data  reduction  program 
that  performs  the  above  experiment  is  described  in  this 
section. 

Figure  22  depicts  the  flow  chart  for  the  Solid  Fuel  Ram- 
jet 'SFRJ'  option,  which  is  assigned  to  special  function  key 
#3  on  the  HP-85.   This  option  is  executed  after  having  loaded 
and  run  the  'HEATER'  program  (figure  14),  completed  the 
'SETUP'  option,  and  executed  the  'FLOW  option  to  set  up  the 
proper  flow  rates  of  air,  oxygen,  and  ethylene.   If  the 
vitiated  air  heater  is  to  be  used,  the  'START'  option  (fig- 
ure 16)  is  executed  to  bring  the  heater  up  to  stable  operat- 
ing temperature.   Once  the  air  temperature  and  total  air  flow 
rate  (m,)  are  within  desired  limits  of  the  required  values, 
the  heater  maintenance  ' HTRMNT '  (figure  18)  is  automatically 
set  to  interrupt  the  computer  at  10  second  intervals  to 
check  the  operation  of  the  heater. 

Once  the  'SFRJ'  key  is  pressed,  the  operation  of  the 
Solid  Fuel  Ramjet  is  initiated.   Both  air  flow  control  valves 
are  initially  opened  to  approximately  the  50%  position. 


66 


09 


--09 


a3 

u 

a 

c 

1 — 

zv 

OJ 

> 

't-^ 

1 — 

o^ 

• 
CO 

> 

( — 1 

2=: 

CD 

OJ 

CO 
03 

U-i 

0R^ 

Ol 

L. 

U-l 

Ol- 

t       ) 

J — 

o 

—  02 


4^ 


—  01 


0 


CD 


C9 


ca 

OJ 


r^a 


CO 


(°(D  a^nssBHd  ffiaHVHO 


67 


AIR-OFF 
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Figure    22.       Solid   Fuel    Ramjet    Flow  Chart 
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WAIT 
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T 


1)  IGN — OFF 

2)  C2H4— OFF 

3)  Pj, — Meas  DVM 


ON 
Tli'IER  in      ] 


▼  YES 

1) 

Meas  ?ti  '  ''•ti  comp  m 

2) 

Meas  P|-2>  T^^  comp  m  ,  m. 

3) 

Meas  Th,  Ph,  T^^,  T^^^ 

4) 

T  =  (mp.Th+ih^p.T^J/ni^ 

5) 

Store  time,  m^,  T,  T^^^ 

SFRJ  Shutdown 

1)  Air — OFF 

2)  N2— ON  (5  sec) 

3)  Heater  Shut 
down 


Figure  22  (Continued) 
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Initially,  heater  output  air  is  routed  via  primary  and  bypass 
air  lines  to  air  dump  valves,  which  are  held  in  the  dump  po- 
sition, venting  heater  air  to  atmosphere.   The  SFRJ  igniter 
and  the  ethylene  valves  are  also  maintained  off  at  this  time. 

The  primary  valve  position  is  then  adjusted  by  the  com- 
puter until  the  primary  air  flow  rate  (m  )  is  within  5%  of 
the  desired  value.   The  computer  programs  the  6  942A,  TIMER/ 
PACER,  A/D  converter,  and  memory  cards  to  read  the  SFRJ  motor 
chamber  pressure  (P-)  transducer  output  voltage  continuously 
for  60  seconds  after  being  externally  triggered.   The  reading 
rate  is  fixed  at  approximately  16  readings  per  second. 

When  the  6  942  is  programmed  the  SFRJ  manual  ignition 
switch  is  enabled  by  the  computer  which  then  monitors  this 
switch  until  it  detects  closure.   Once  the  switch  is  closed 
and  sensed,  the  Timer /Pacer  card  is  triggered  and  the  com- 
puter timer  is  set  to  zero.   Heater  air  is  then  directed 
into  the  SFRJ  motor,  the  igniter  is  turned  on,  and  the  C^H 
is  turned  on. 

The  computer  waits  a  preset  time  in  seconds  (T    )  before 
shutting  off  the  igniter  circuit  and  the  ^2^1^   valve.   During 
the  start  sequence  the  3456A  digital  voltmeter  is  set  to 
monitor  chamber  pressure  (Pp) .   The  chamber  pressure  is  mea- 
sured and  compared  to  the  desired  steady  state  value.   If 
this  pressure  exceeds  .75  of  the  desired  value  by  the  end  of 
ignition  the  computer  proceeds  to  a  measurement  sequence  of 
several  pressures,  temperatures,  and  flow  rates  until  the 
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desired  burn  time  in  seconds  expires.   However,  if  the  pres- 
sure is  less  than  .75  of  the  desired  value  the  program  di- 
rects the  user  back  to  attempt  another  ignition  sequence. 
The  automatic  measurement  sequence  proceeds  from  the 

point  where  the  pressure  exceeds  .75  P^  ■,   .   ^  as  (refer 
^  '^  C  desired 

to  figure  21) : 

a)  Pt-i  /  Tt-,  are  measured 

b)  m  .   =  

air        Q 

where  C_^  =  .97 

A  =  TTd^/4 
d  =  .235 


'.  =j'(*P  * 


Y 
Y  =  1.4 


P-Pp  /  Trp   measured  at  primary  air  orifice 
PT9  ^9  C, 


d)   mp  = 


[R  Tt2 
where  Q  =/ 

e)   V™   is  set  to  obtain  desired  m 
n^p  p 

•  •         • 

^^   "^BYPASS  ^  ^A  ~  ^PRIMARY 

^BYPASS  =  0  if  ^o  bypass  selected 


^^   ^HOT  AIR  ^e^sured 
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h)   T   measured 

i)   SFRJ  head  pressure  and  temperature  measured  P„,  T„ 

n    n 

^'  ^  =  s- 

air 

k)   Values  of  time,  mair^  ^'  ^^^  ^hot  ^^^  stored  in  an 
array  for  later  use  in  reduction  of  data. 

The  above  measurements  proceed  for  a  total  of  100  read- 
ings at  a  rate  of  approximately  1.5  per  second. 

Both  the  SFRJ  and  the  heater  are  shut  down  automatically 
when  timer  #1  expires  at  (t  )  seconds  after  successful 
ignition. 

Once  the  measurement  sequence  completes  and  the  shutdown 
sequence  is  completed  the  computer  then  waits  until  the  multi- 
programmer  memory  board  signals  completion  of  its  storage 
operation. 

During  the  operation  of  the  SFRJ  option,  special  function 
key  #4  is  assigned  to  the  SFRJ  shutdown  routine.   This  pro- 
gram is  automatically  executed  t,  seconds  after  ignition, 
however,  it  can  be  invoked  at  any  time  by  the  operator  and 
when  called  performs  the  following  steps: 

a)  Air  valve  to  dump  position 

b)  Nitrogen  purge  on  for  5  seconds 

c)  Heater  shutdown  sequence. 

After  all  timed  events  have  completed  the  computer  auto- 
matically loads  and  executes   'SFRJ2',  which  is  the  data 
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reduction  program  for  the  Solid  Fuel  Ramjet  experiment.   Fig- 
ure 23  shows  the  flowchart  of  the  data  reduction  program. 

A  large  quantity  of  measured  data  is  automatically  passed 
from  the  Heater  program  to  the  reduction  program.   The  com- 
puter then  prompts  the  user  for  other  required  input  data, 
such  as: 

a)  fuel  type   HTPB/PMM 

b)  final  grain  weight  (grams) 

c)  initial  grain  weight  (grams) 

d)  initial  grain  internal  diameter  (cm) 

e)  final  aft-end  internal  diameter  (cm) 

f)  grain  length  (cm) 

2 

g)  area  of  the  exhaust  nozzle  throat  (cm  ) 

With  this  data  entered  and  stored,  the  computer  then  reads 
in  the  chamber  pressure  (Pp)  array,  which  is  a  1000  element 
vector  of  transducer  voltages.   These  voltages  are  subsequent- 
ly converted  to  pressure  (PSIA)  by  the  formula 

P^(I)  =  K^(V(I)  -  V0)  +  B9 

where    K^  =  79.4  34 

V(I)  =  transducer  voltage 

V0  =  initial  transducer  voltage  at  zero  psig 
B9  =  atmospheric  pressure 
The  computer  then  locates  the  maximum  pressure  reading 
and  its  location  in  the  array.   The  computer  then  searches 
both  forward  and  backward  in  time  from  the  location  of  the 
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INPUT  DATA 

1)  Enter  Init.  Grain  Wt. 

2)  Enter  Final  Grain  Wt. 

3)  Enter  Grain  Length 

4)  Enter  Grain  Diameter 

5)  Enter  Final  Aft  Diam. 

6)  Enter  Fuel  Type 


1)  Read  in  P(,(t) 

2)  Convert  P(t)  to  PSIG 


1)  Find  Pp(inax) 

2)  Find  Tine  (Pc  max) 

3)  FindTs=Tiine(.2.p^  ^^) 

4)  FindTv=Tine(.2^c  max) 

5)  Tk  =  Tt7  -  Tc 


^C  = 


INTEGRATE 
Tp 

/      P    (t)dt 
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Figure    23.      Solid  Fuel   Ramjet   Reduction   Flow  Chart 
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Figure    23     (Continued) 
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maximum  pressure  until  it  finds  the  time  where  pressure  has 

dropped  to  0.2  times  the  maximum  pressure  and  records  these 

locations  as  t  ^   ^  and  t^.  .  ,  .   The  effective  final  burn 
start      finish 

time  is  t_,_  =  t_  -  t^ .   The  computer  then  performs  a  numerical 

Br      r      b 

integration  of  the  pressure  time  trace  from  t^  to  t„  and 
divides  the  result  by  t„p  to  give  Pp ,  the  average  pressure. 

/  Pp(t)dt 
t 

^c  =  -^— E 

Now  that  t   and  tp  are  known,  the  computer  can  locate 

these  times  in  the  time  column  of  the  data  array  (100x3 

matrix)  and  note  the  row  numbers  of  I   and  I  .   The  average 

air  temperature  T    and  flow  rate  lii  .   are  then  computed 
^         ave  air  ^ 

between  t<^  and  t„   by: 

I, 


F 

m 

A 


L* 


rti 


air 

(ip- 

^s' 

E 

T 

i=is 

"^av    (1^,-13) 


Data  reduction  can  now  proceed  as  follows 


AW.  =  Initial  W   -  Final  W 


4  AW 
d^  =   f  +    d. 
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3 
where  pp  =  Fuel  propellant  density  (gms/cm  ) 

L  =  length  .  (cm) 

d.  =  initial  diameter  (cm) . 

The  average  burn  rate  is  then 


d.  -  d. 
f     1 

r     =  -rr- 

ave     2tgp 


The  equivalent  fuel  flow  rate 

AW. 
t 

BF 
Thus,  fuel/air  and  equivalence  ratios  are 


and 


(j)  = 


^STOICH 


Then  with  T,  't> ,    and  the  type  of  fuel  burned,  the  theoretical 
values  of  T4-.^^/  Rp ,  and  v^  are  computed  from  a  polynomial 
regression  obtained  from  the  PEPCODE  [Ref.  5]  computer  pro- 
gram at  the  Naval  Postgraduate  School  computer  center.   The 
effect  of  P_  on  T-j-,  u  is  negligible  for  nominal  variations 
(40-60psia) .   However,  Pp  effects  can  be  readily  incorporated 
if  desired.   These  polynomial  equations  are  given  in  Appendix 
B. 
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Once  the  values  of  T-|-,r^,  Rp  /  and  y   are  known  then 


and  the  experimental  value  of  chamber  stagnation  temperature 

T-i-     is  determined  from 
^exp 

X  2 

^      _  .  ^C  ^T  ^D  ^2 

^exp   .  , .     ,  -  , 
"^^   \  m    +  m„ 
\   ave    F 

where  A^  =  throat  area 
C   =   97 
Finally,  the  efficiency  is  computed 

T*.  -  T 


-exp 
^tth 


-  T 


Hard  copy  printouts  of  all  data  is  then  provided  by  the 
system  printer. 
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APPENDIX  A 
LIST  OF  HP-IB  ADDRESSES 

1  CRT  (internal) 

2  Printer  (internal) 
701  2631B  Line  Printer 
705  7225  Plotter 

709  3497A  Data  Acquisition/Control  Unit 

721  HP-85  Controller 

722  3456A  High  Accuracy  Digital  Voltmeter 

723  6942A  Multiprograimner 

724  3437A  High  Speed  Digital  Voltmeter 
HG-78  Scanivalve  Interface  #1 
HG-78  Scanivalve  Interface  #2 

":D700"  Disc  Drive  #0 
":D701"  Disc  Drive  #1 
":T"  Internal  Tape  Drive 
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APPENDIX  B 
PEPCODE  REGRESSION  EQUATIONS  FOR  HTPB  AND  PMM 


HTPB  =  Fuel  Type        ^STOICH  "  '^"^^        '^P  ^  -0332  #m/in^ 


(})  =  := P-  =  40-60  psia 

STOICH 


Tt    =  -933.3  +  4893^  +  203(j3^  -  Q99<t>^    -  205.90^  +  1.1925T^ 


+    1.014(10"'^)T^^  -  1.84(10  ^)T^^  -  1.138T^   +  .3864T^^^ 
-  7.67(10"^)f^^(j) 


Y^  =  1.40162  -  .277234)  +  .14834>^  -  .02448(1)^  +  .005234)^ 

-  4.994(10"^)T-  +  4.255(10"^)T^^  -  2  .  36  ( lO""^^)  T^  ^ 

A  A  A 

+  7.7585(10"^)T.(|)  -  3.1138  (10"^)f.'t)^  -  1 .  885  ( lO"^)  f^^  4) 

A  Pi  A 


R^  =  52.966  +  3.516<J)  -  9.55^)^  +  5.71(|)^  +  .2739(|)^ 

-  S.OSdO""^)!^^  +  1.36(10"^)T^^  -  5.99(10"-^^)T^^ 

+  1.0208(10"^)T^(j)  -  1.872  (10"^)T^(t)^  +  2  .  185  ( lO""^)  T^^4) 
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PMM  =  Fuel  Type        ^STOICH  =  ^'^^        pp  =  .0426  #m/in^ 

P   =  40-60  psia 


T^         =    1642  -  5902^  +  14688^^  -  7478(|)^  +  4.05T^(]3 

-  .18f^  -  4.68T^(j)^  +  1.28T^4>'^  -  3 .  259  ( lO"^)  T^^(|)^ 


Y^  =  1.34  -  .0265(|)  -  9.126  (10"^)T^  -  .193(t)^  +  .1329(t>^ 

-    4.205  ( 10"^)  f^(|)  +  7.429  (10"^)f^(|)^  -  2  .12  {10~^)T^(t>^ 
+  2.09  (10~^)T^'^({)^ 


R^=  47.1569  +  ^[2  7.54  +  ^{2.9108-(f)(41.242-(})[l9.54  5-11.56  7(J)^]))J 
+  (10"'^)f^|l2403  -  3.2029  (10"^)T^^ 
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APPENDIX  C 
PROGRAM  LISTINGS 


iO  •  HELP  FOR  HEATER 
20  OPTION  BASE  i 

30  COM  y(iOO,4) ,M9,N0$[i] ,P0,V9,B9 
40  CRT  IS  2 

SO  DISP  "The  HEATER  prograw  contains  project  control  an 
d" 

60  DISP  "  DATA  ACQUISITION  and  DATA  REDUCTION  prograMS 
for  the  VITIATED  AIR  HEATER." 

70  DISP  "PRESS  'SET  UP'  KEY  TO  MAKE  PRELIMINARY  AIR  FLO 
U  CALCULATIONS. " 

80  DISP  "WHEN  YOU  ARE  SATISFIED  WITH  YOUR  CALCULATIONS 
YOU  CAN  FLOy  THE  AIR  SYSTEM  TO  CHECK .  " 
90  DISP 

iOO  DISP  "PRESS  'FLOy  KEY  TO  COLD  FLOW  THE  AIR  TO  MEAS 
URE  Pti  AND  Tti  AND  COMPUTE  MidotiAIR." 

iiO  DISP  "IF  YOU  ARE  NOT  SATISFIED  WITH  Pti  SETTING  OR 
ACTUAL  AIR  FLOW  THEN  'FLOW'  IS  RUN  AGAIN." 
120  DISP  "AFTER  'FLOW'  IS  EXECUTED  THE  'START'  OPTION  K 
EY  IS  PRESSED  TO  START  HEATER  SEQUENCE." 

130  DISP  "'START'  BRINGS  HEATER  UP  TO  DESIRED  TEMPERATU 
RE  .  " 

140  DISP  "WHEN  DESIRED  TEMP  IS  ATTAINED  THE  'HTRMNT'  HE 
ATER  MAINT  ROUTINE  IS  RUN  EVERY  20  SECS . " 
ISO  DISP  "THE  HEATER  IGNITER  IS  DE-ENERGIZED  AFTER  8  SE 
CONDS  OF  OPERATION  OR  WHEN  700  DEG  IS  " 

160  DISP  "REACHED  ,  ALSO  IF  THE  AIR  TEMP  DROPS  BY  50  DE 
G  IN  5  SECONDS  DURING  START  UP  THE  IGN" 
170  DISP  "  IS  SHUT  OFF  AND  THE  HEATER  IS  SECURED." 
180  DISP  "WHEN  THE  AIR  TEMPERATURE  IS  WITHIN  ISO  DEGREE 
S   OF  THAT  DESIRED,  THE  HEATER  MAINT  " 

190  DISP  "  ROUTINE  IS  CALLED.  THIS   PROGRAM  EXECUTES  E^ 
ERY  20  SECONDS  UNTIL  ANOTHER  OPTION" 
200  DISP  "  IS  SELECTED,  MEASURING  AND  RECORDING  THE  AIR 

TEMP  .  " 
210  DISP  "  THE  SOLID  FUEL  RAMJET  OPTION  KEY  'SFRJ'  IS  P 
RESSED  TO  ENERGIZE  THE  SFRJ  WITH  OR  " 

220  DISP  "WITHOUT  THE  HEATER  IN  OPERATION.  THIS  OPTION 
RUNS  FOR  A  PRESET  TIME  BEFORE  SECURING" 

230  DISP  "  'SFRJ2'  THE  DATA  REDUCTION  PROGRAM  IS  AUTOMA 
TICALLY  CALLED  AND  EXECUTED  AFTER  'SFRJ'" 
240  CRT  IS  1 
250  CHAIN  "HEATER" 
260  END 
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iO  !  HEATER 

20  OPTION  &ASE  i 

30  COM  V(iOO,4) ,M9,N0$[tl ^PO,^?,^? 

40  GN  KEY*  i,"SET  UP"  GOSUB  i'SQ 

50  ON  KEY*  2,"FL0U"  GOSUB  790 

60  ON  KEY*  3,"SFRJ"  GOTO  2270 

70  ON  KEY*  7,"TJTC"  GOSUB  1220 

80  ON  KEY*  S,"HELP"  GOSUB  1240 

90  ON  KEY*  6/'START"  GOSUB  1270 

100  ON  KEY*  4, "SHUTDOWN"  GOSUB  2140 

110  CLEAR  e  KEY  LABEL 

120  DISP  "SELECT  OPTION" 

130  GOTO  130 

140  !  Initialize  the  constants  as  needed  to  do  the  prel 

irtinary 

150  !  calculations  for  setting   up  the  air  heater. 

160  DISP  "UILL  YOU  BE  RUNNING  THE  HEATER  ONLY?  <Y/N)' 

170  INPUT  Y$ 

180  IF  Y$="N"  THEN  GOTO  2880 

190  ! 

200  OUTPUT  709  j "D04 , 4 ,5 , 6 >7" 

210  LET  ^^9=280  !  STORAGE  TANK  VOLUME 

220  G=1.4  !  AIR  GAMMA 

230  DISP  "ENTER  STORAGE  TANK  GAGE  PRESSURE  (PSIG)" 

240  INPUT  B8 

250  LET  C=.97  !  DISCHARGE  COEF . 

260  LET  Gl  =  32.2  !  GRA^^ITY 

270  LET  D=1.5  !  INCHES 

280  Dl=.235  !  Dstar  choked  ventur.L 

290  LET  R=S3.3  !  GAS  CONSTANT  FOR  AIR 

300  CLEAR 

310  DISP  "ENTER  DESIRED  MASS  FLOW  RATE  (*/SEC)  M9 " 

320  INPUT  M9 

330  PRINT  "  FOR  DESIRED  MASS  FLOW  RATE  OF  ";M9;"  */SEC . 
It 

340  !  DISP  "ENTER  NOZZLE  DIAMETER    (Di=COLD  AIR=D*)" 

350  !  INPUT  Dl 

360  DISP  "ENTER  COLD  AIR  TEMPERATURE  <DEG==F)  Tl" 

370  INPUT  Ti 

380  LET  Ti=Ti+460 

390  DISP  "ENTER  DESIRED  HOT  AIR  TEMPERATURE  (DEG=R)  T8" 

400  INPUT  T8 

410  DISP  "ENTER  BAROMETRIC  PRESSURE  SETTING  (INCHES  HG ) 

(29.92«STANDARD  SEA  LEVEL)" 
420  INPUT  B9 
430  B9=.4913*E<S' 
440  !  B9  WILL  BE  USED  TO  CONVERT  ABS  PRESSURES  TO  GAGE" 

450  !  NOW  COMPUTE  VISCOSITY  FOR  LATER  USE" 
460  LET  U=.0000000229*T8*1 .5/(T8+198.6) 
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470 

1 

480 

!  COHPUTE  FUEL  AIR  RATIO  FROM  PEF 

490 

LET  F=<T8-550)/78000 

500 

t  , 

5i0 

!  This  F  was  for  air  at  approx.  ' 

!  co«pute  ethylene  flow  rate 

520 

LET  M8=F*M9  !  MDOT  FUEL 

530 

M7=96/28*M8  !  OXYGEN  FLOW  RATE 

540 

V=i83720*h7  !  OXYGEN  SET  VOLTAGE 

550 

V$=VAL$(V) 

560 

OUTPUT  709  ;"AO3,0/'^V1 

570 

CLEAR 

580 

DISP  USING  590  ;  M8;826.06*M8 

590 

IMAGE  -SET  ETHYLENE  FLOW  RATE=-, 

D.3D."  SCFM" 

)DE  EXCERT 
70  degrees  farenhei 


ID.SD,"  */SEC">5X,;;> 

600  DISP  USING  610  j  M7;722.58*M7 

610  IMAGE  "SET  OXYGEN  FLOW  RATE=" , ID . 5D /'  */SEC " , 5X, 2D. 

3D/'  SCFM" 

620  PRINT  "OXYGEN  VOLTAGE=  "jV;"  MILLIVOLTS" 

630  LET  A=PI/4*Di^2 

640  LET  F1=SQR<G*<2/<G+1) )^( (G+l)/<G-i) ) ) 

650  PRINT  "F  < GAMMA ):=";Fi 

660  LET  Q9=SQR(Ti*R/Gl) 

670  LET  P6=M9*Q9/(C*A*Fi)  !  PSIA 

680  LET  P5==P6-B9  •  PSIG 

690  DISP  USING  70  0  \    P  V. 

700  IMAGE  "SET  PT  TO  ",6D/'   PSIG" 

710  LET  T5=<B8-P5*2>*V9*144/<R*T1*M9) 

720  DISP  USING  730  ;    T5 

730  IMAGE  "RUN  TIME  IS  ",5D.D,-  SECONDS' 

740  COPY 

750  CLEAR 

760  DISP  "SELECT  OPTION" 

770  KEY  LABEL 

780  RETURN 

790  !  FLOW  OPTION 

800  OUTPUT  709  ;  "DQ4  ,4 ,5,6 , 7  "  !  AIR  VALVE  CLOSED 

810  CLEAR 

820  WAIT  2000 

830  !  THIS  SEGMENT  FLOWS  THE  SYSTEM  COLD  TO  DETERMINE  A 

CTUAL  MASS  FLOW  RATE 

840  OUTPUT  709  ;"AC5"  !  MEASURE  PTl 

850  OUTPUT  722  ; "HSM002SW2Z0SO1L1 . ISTISOF 1R1T3QX1  ' 

860  GOSUB  2230 

870  VO=V 

880  PRINT  "VO=",VO 

890  OUTPUT  709  ;"DC4,6"  !  AIR  VALVE  OPENS 

900  WAIT  iOOO 

910  OUTPUT  722  ;  "HSM002SW2Z0SO1L1 . ISTISOFIR iT3QXl " 

920  GOSUB  2230 
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930  PRINT  "^=",y 

940  P=79.434*(y-V0)  !  VOLTS  TO  PSIG 

9S0  PRINT  "PTi=",P 

960  OUTPUT  709  }"AC6"  !  MEASURE  TTi 

970  OUTPUT  722  ; "HSW2SM00220SOlLi . iSTISOFiR ITiQX i  " 

980  GOSUB  2230 

990  T=4i57S*y+S36.6  !  VOLTS  TO  DEGREES  R 

iOOO  OUTPUT  709  j  "004,6"  !  CLOSE   AIR  VAL^E: 

iOiO  !  CONVERSIONS  ON  T  AND  P  NEED  TO  BE  DONE. 

1020  Q9=SQR<R*T/Gi) 

1030  M=C*A*F1*(P+B9)/Q9  !  ACTUAL  MA  FLOW 

1040  DISP  "ACTUAL  MASS  FLOW  RATEE  IS    =";M;"  *M/3ec  .  " 

1050  M8=F*M  e  M7=96/28*M8 

1060  V=i83720*M7  !  OXYGEN  SET  VOLTAGE 

1070  OUTPUT  709  j "A03 , 0 /'iVAL* ( V ) 

1080  P7=M9*Q9/<Fi*C*A)  !  PSIA 

1090  P=P7-P-B9  !  PRESSURE  INCREMENT 

1100  PRINT  "  DELTA  PTi=",P 

1110  DISP  USING  1120  \    P,M9 

1120  IMAGE  "CHANGE  Ptl  BY  ",4D,"  PSIG  TO  ACHEIVE  A  Mdot 

(AIR)  OF  ",D.5D,"  *«ass/sec." 
1130  DISP  "ARE  YOU  SATISFIED  WITH   THE  MASS  FLOW  RATE  A 
ND  DOME  PRESSURE  SETTING?  (Y/N)" 
1140  INPUT  Y* 

1150  IF  Y$="Y"  THEN  GOTO  1180 
1160  T1=:T  !  P6=P7!  SAVE  LAST  PTl  AND  TTi 
1170  GOTO  790 
1180  DISP  "SELECT  OPTION' 
1190  M9=M  !  ACTUAL=DESIRED  MDOT 
1200  KEY  LABEL 
1210  RETURN 
1220  BEEP  200,100 
1230  RETURN 
1240  !  HELP 
1250  CHAIN  "HELP" 
1260  RETURN 

1270  !  HEATER  START  UP  ROUTINE 
1280  OUTPUT  709  ;  "004,6"  !  AIR  ON 
1290  CLEAR 

1300  OUTPUT  709  ;"AC5"  !  Ptl 

1310  OUTPUT  722  ; "HSW2SM00220SQ1L1 . ISTISOFlRlTlQXi" 
1320  GOSUB  2230 

1330  P=79.434*(V-V0)  !  VOLTS  TO  PSIG 
1340  OUTPUT  709  ;"AC6"  !  Ttl 

1350  OUTPUT  722  ; "HSW2SM002Z0SO1L1 . ISTISOFIRITIQXI " 
1360  GOSUB  2230 
1370  T=41575*V+536.6  !  T  COLD 
1380  Q9=SQR(R*T/G1) 
1390  M=C*A*F1*(P+B9)/Q9  !  ACTUAL  MDOT  AIR 
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1400  M8=F*h  9    M7=96/28*h8 

1410  ^=i8372Q*h7  !  02  SET  ^OLTS 

1420  OUTPUT  709  ;  "A03  ,  0 /'A,yAL$  <  U  ) 

1430  OUTPUT  709  ;"AC7"  !  T  hot 

1440  OUTPUT  709  ;"DC4,7"  !  IGNITER  OM 

14S0  OUTPUT  709  ;"DC4,4,5"  !  02  iC2H4  ON 

1460  ON  TIMER*  1,8000  GOSUB  1670 

1470  T3=T  !  T  LAST 

1480  WAIT  5000 

1490  OUTPUT  722  j "HSU2SM002Z0SOlLi . ISTISOFIR ITIQXI " 

1500  GOSUB  2230 

1510  T2=i002.3*V+457.53  !  Thot 

1520  Tl=T2-T3+50 

1530  IF  T2>700  THEN  GOSUB  1720 

1540  IF  TKO  THEN  1760 

1550  PRINT  "T!hot=",T2 

1560  I  CHECK  TO  SEE  IF  HEATER  IS  WITHIN  DESIRED  LIMITS 

1570  IF  T2+150-T8>0  THEN  1840 

1580  !  HEATER  IS  NOT  YET  NEAR  OPERATING  RANGE 

1590  OUTPUT  709  ;  '•AC7" 

1600  T3=T2  !  LAST  T  HOT 

1610  WAIT  4000 

1620  OUTPUT  722  ; "HSW2SM002Z0SO1L1 . ISTISOFIR ITlQXi  " 

1630  GOSUB  2230 

1640  T2=1002.3*V+457.53  !  Thot 

1650  Tl=T2-T3+50 

1660  GOTO  1540  !  READ  T  HOT  AGAIN 

1670  !  SHUT  OFF  IGNITER 

1680  OUTPUT  709  j"D04,7" 

1690  PRINT  "IGNITER  SHUT  DOWN-5  SECS . " 

1700  OFF  TIMER*  1 

1710  RETURN 

1720  OUTPUT  709  }"D04,7" 

1730  OFF  TIMER*  1 

1740  DISP  "  T:HOT  >  700  R  IGNITER  OFF" 

1750  RETURN 

1760  !  TEMP  GRAD  IS  NEG  TIME  FOR  SHUT  DOWN 

1770  OUTPUT  709  ;"D04,7" 

1780  OUTPUT  709  ;"D04,4,5" 

1790  OUTPUT  709  ;"AO3,0,0" 

1800  OUTPUT  709  ;"D04,6"  !  AIR  OFF 

1810  DISP  "  HEATER  SHUTDOWN-NEG  TEMP  GRADIENT* 

1820  OFF  TIMER*  1 

1830  STOP 

1840  I  HEATER  IS  NEAR  OPERATING  POINT 

1850  OFF  TIMER*  1 

1860  PRINT  "  HEY  I'M  CLOSE" 

1870  ON  TIMER*  3,10000  GOSUB  1940 

1880  DISP  "ENTER  HEATER  MAINT  ROUTINE" 

1890  CLEAR 
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1900  DISP  "HEATER  RUNNING" 

1910  UAIT  500 

1920  KEY  LABEL 

1930  GOTO  1890 

1940  !  HEATER  MAINTENANCE 

1950  OUTPUT  709  ;"AC7" 

1960  OUTPUT  722  ; "HSM002SU2Z0SO1L1FL0  .  0 ISTISOFIRITIQXI  " 

1970  GOSUB  2230  !  TAKE  READING 

1980  T2=1002.3*U+457.53  !  Thot 

1990  PRINT  '•T:hot:AIR="jT2 

2000  T4=T8-T2  !  DELTA  T:hot 

2010  IF  ABS(T4)>200  THEN  GOTO  2140 

2020  IF  ABS(T4)<100  THEN  1940 

2030  IF  T4<0  THEN  GOTO  2060 

2040  DISP  ••  NOT  HOT  ENOUGH  " 

2050  GOTO  2080 

2060  DISP  "  TOO  HOT  " 

2070  GOTO  2080 

2080  F4=T4/7a000 

2090  F=F+F4  !  NEW  FUEL  AIR  RATIO 

2100  M8=F*M9  8  M7=96/28*M8 

2110  V=183720*M7 

2120  OUTPUT  709  ;  "A03 , 0 /'^yAL$ ( V) 

2130  RETURN 

2140  !  HEATER  SHUTDOWN 

2150  OUTPUT  709  }"AO3,0,0' 

2160  OUTPUT  709  ;"D04,7,4,5" 

2170  PRINT  "HEATER  SHUTDOWN" 

2180  PRINT  "  LAST  T:hot=";T2 

2190  WAIT  5000 

2200  OUTPUT  709  j "D04 , 4 ,5 ,6 ,7" 

2210  OFF  TIMER*  3 

2220  RETURN 

2230  !  SRQ  SERVICE  ROUTINE 

2240  STATUS  7,1  ;  Z9 

2250  ENTER  722  j  V 

2260  ENABLE  INTR  7;8  0  CLEAR  722  8  RETURN 

2270  !  SOLID  FUEL  RArtJET  OPTION 

2280  OUTPUT  709  ; "D04 ,8 , 9 , 10 , i 1 , 12 , 13"  !  AIR  DUMP 

2290  OFF  KEY*  4  8  ON  KEY*  4,  "SHUTDOWN"  GOSUB  3050 

2300  IF  N0$="N"  THEN  2330 

2310  OUTPUT  709  ;"DC4,9,13" 

2320  DIM  O0$C90] 

2330  O0$="WF,3. 1,0/3.2,0,3.3,0T,SF,2,3,1,. 00 1 , 12T ,MI  ,2 , 

1 0  0  OT , WF , 3 , 999T , AC , 3T , WF , 4 . 2 , 1 T " 

2340  OUTPUT  723  ;O0$ 

2350  OUTPUT  723  ' "OP  4  . 03ST" 

2360  CLEAR  722  8 'oUTPuf  722  ; "HSM002SU2Z0SO1L1FL0 . OlSTI 

S0FiRlT3Q" 
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2370  OUTPUT  709  ;"DC4,i4"  ©  DISP  "READY  FOR  SFRJ  IGNITI 

ON" 

2380  OUTPUT  709  ;"ACi3"  !  Pchanber 

2390  OUTPUT  722  j"Xl"  »  GOSUB  2230 

2400  y9=g  !  INITIAL  Pc  READING 

2410  OUTPUT  709  ;"ACi4"  !  IGNITION 

2420  OUTPUT  722  ;"Xi"  8  GOSUB  2230 

2430  IF  ABS(V)<3  THEN  2420 

2440  !  BEGIN  SFRJ  START  UP 

2450  SETTIHE  0,0 

2460  OUTPUT  709  ) "DC4 , 8 , i i , 12" 

2470  WAIT  i000*T6 

2480  OUTPUT  709  j"D04,i2,ll" 

2490  OUTPUT  709  )  "AC13"  !  Pchanber 

2500  OUTPUT  722  ;"X1"  S  GOSUB  2230 

2510  IF  y> .75»(P0/151.5)  THEN  2540 

2520  CLEAR 

2530  DISP  "NO  IGNITION  TRY  AGAIN"  S  GOTO  2420 

2540  !  BEGIN  MEASUREMENT  SEQ. 

2550  ON  TIMER*  i,T6*1000  GOSUB  3050 

2560  OUTPUT  722  ; "HSM002Sy2Z0SQlLlFL0 . OiSTISOFiR iT3«" 

2570  FOR  J=l  TO  iOO 

2580  OUTPUT  709  ;"AC5"  e  OUTPUT  722  ;"Xi"  9    GOSUB  2230 

2590  P=79.434*(V-V0)  8  OUTPUT  709  i"AC6"  9    OUTPUT  722  j 

"XI"  8  GOSUB  2230 

2600  T  =  41575*V+536.6  8  ^;(  J  ,  1  )=C*A*Fi*<P+B9) /SQR  (R*T/G1 ) 

2610  OUTPUT  709  i"ACi2"  8  OUTPUT  722  ;"Xi"  S  GOSUB  2231 

2620  T3=41575*g+536.6  8  OUTPUT  709  j"AC7"  8  OUTPUT  722 
;  "XI"  8  GOSUB  2230 

2630  V(J,4)=i002. 3*^+457. 5  8  T2=V<J,4)  !  T:hot 
2640  OUTPUT  709  ; "AC8"  8  OUTPUT  722  ;"Xi"  8  GOSUB  2230 
2650  P-79.434*(y-gO)  8  OUTPUT  709  j"AC9"  8  OUTPUT  722  ; 
"Xi"  8  GOSUB  2230 

2660  T=4i575*V+536.6  8  M3=C*A*F1*<P+B9) /SQR <R*T/Gi ) 
2670  M4=V(J,i)-M3  8  <J(J,3)=TIME 
2680  IF  N0*="N"  THEN  M4=0 
2690  OUTPUT  709  ;"ACil"  8  OUTPUT  722  ;"Xi"  8  GOSUB  2230 

2700  P 4=79.  434* (V-yO)  8  OUTPUT  709  ;"AC10''  8  OUTPUT  722: 

;"Xi"  8  GOSUB  2230 

2710  T4=41575*y+536.6  8  ^ ( J ,2 )=(M3*T4+M4*T3)/V< J , 1 )  !  T 
bar 

2720  PRINT  "J=";J,  "TIME=:")V<  J,3) 

2730  PRINT  "Thead=",T4 

2740  PRINT  "Pheacl  =  ",P4 

2750  PRINT  "T:hot  =  ",'v'(J,4) 

2760  PRINT  "Midot  AIR=",y<J,i) 
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2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

Y/N) 

2890 

290  0 

2910 

2920 

T  TH 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 


PRINT  "«:dot  PRI=-",M3 

PRINT  "T;bypass=",T3 

PRINT  "T:BAR=",y< J,2) 

NEXT  J 

O0=SPOLL<723)  !  WAIT  FOR  MEh  INTR 

IF  00<  >64  THEN  2810 

DISP  "hULTIPROGRAHMER  HAS  P:chartber  data  " 

SEND  7  ;  UNL  MLA  TALK  23  SCG  12 

ENTER  723  ;  Oie  DISP  "  SLOT#="j01 

OUTPUT  723  ; '•DC>3T,UF,4.2,0T,OP,4,0T" 

CHAIN  "SFRJ2" 

DISP  "WILL  YOU  BE  RUNNING  THE  SOLID  FUEL  RAMJET? 


THEN  GOTO  2940 


COMBUSTOR  OPTION  IS  NOT  AVAILABLE  A 
ANOTHER  OPTION" 


PRESSURE  (PSIA)" 
TIhE   IN  (SECONDS)" 


INPUT  Y* 
IF  Y$="Y" 
CLEAR 

DISP  "TURBOJET 
IS  TIME.  SELECT 
GOTO  150 
!  SFRJ  OPTIONSi 
CLEAR 

DISP  "ENTER  DESIRED  CHAMBER 
INPUT  PO 

DISP  "ENTER  THE  DESIRED  BURN 
INPUT  T5 

DISP  "ENTER  THE  IGNITION  TIME 
INPUT  T6 

DISP  "WILL  THE  BYPASS  OPTION  BE  USED? 
INPUT  N0$ 
GOTO  190 
!  SFRJ  SHUTDOWN 

OUTPUT  709  }"D04,8>9,11,12,13" 
OUTPUT  709  ;"DC4,10"  !  N2  PURGE  ON 
WAIT  5000  e  OUTPUT  709  ;"DO4,10" 
GOSUB  2140  !  HEATER  SHUTDOWN 
CLEAR  &    DISP  "SFRJ4  HEATER  SHUTDOWN  " 
RETURN 


IN  (SECONDS)" 
(Y/N)" 


9    KEY  LABEL 
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iO  !  SFRJ2  DATA  REDUCTION 

20  OPTION  BASE  i 

30  COM  U(iOO,4)  ,M9,N0$[il  ,P0,iv'7,B9 

40  SHORT  P(IQOO) 

SO  GCLEAR 

60  CLEAR 

70  Gi^32.174  !  gc 

80  C=.97  !  DISCHARGE  COEF . 

90  DISP  "SFRJ  DATA  REDUCTIONS" 

iOO  DISP  "ENTER  INITIAL  GRAIN  WEIGHT= < GRAMS ) " 

iiO  INPUT  yt 

120  DISP  "ENTER  FINAL  GRAIN  W£IGHT= ( GRAMS) " 

130  INPUT  W2 

140  DISP  "ENTER  GRAIN  LENGTH* (CM ) " 

150  INPUT  L2 

160  L2-^L2/2.54 

170  DISP  "ENTER  INITIAL  GRAIN  DIAMETER=(CM) " 

180  INPUT  DO 

190  D0=D0/2.54 

200  DISP  "ENTER  FINAL  AFT  DIAMETER=(CM) " 

210  INPUT  D2 

220  D2=D2/2.S4 

230  A1=PI*1. 2^2/4  !  AREA  THROAT 

240  DISP  "ENTER  FUEL  TYPE= ( HTPB/PMM ) " 

250  INPUT  Fl$ 

260  IF  F1*="HTPB"  OR  F1$="PMM"  THEN  290 

270  DISP  "ONLY  FUEL  TYPES  AVAILABLE  ARE  <HTPB>  OR  <PMh> 
II 

?.QQ    GOTO  240 

290  U3=<Ui-U2)/453.6  !  DELTA  WEIGHT  IN  *m 

300  I  ENTER  Pc  DATA 

310  SEND  7  ;  UNL  MLA  TALK  23  SCG  5 

320  FOR  1=1  TO  997  STEP  4 

330  ENTER  723  USING  "*,K,K,K,K"  ;  P ( I )  ,P < I  +  i ) , P ( 1+2) , P ( 

1+3) 

340  NEXT  I 

350  ENTER  7 

360  ! 

370  MAT  P=(5)*P 

380  FOR  1=1  TO  iOGO 

390  P(I)=i5.37*P(I)+.22  !  PSIG 

400  NEXT  I 

410  ZO=AMAX(P) 

420  J0=AMAXROU  9  T0=JO*.O6 

430  PRINT  "MAX  Pc=";ZO;"  AT  TIME=";TO;"  SECONDS" 

440  FOR  I=JO  TO  1000 

450  IF  P(I)>.2*Z0  THEN  490 

460  S2=I 

470  PRINT  "S2="iS2 

480  GOTO  500 
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490 

NEXT  I 

500 

FOR  I=JO  TO  1  STEP  -1 

510 

IF  P(I)  >  .2*Z0  THEN  550 

520 

S1  =  I 

530 

PRINT  "S1=";S1 

540 

GOTO  560 

550 

NEXT  I 

560 

S3=S2-S1  !  ELEMENTS  INCLUDED 

570 

34= . 06*S1  e  35= . 06*32 

580 

S6=. 06*33  !  TIME  TBF 

590 

REDIM  P(S3) 

600 

S7=SUM(P) 

6i0 

Pi=S7/33+B9  !  Pc  bar<PSIA) 

620 

REDIM  P(IOOO) 

630 

GCLEAR  e  LOCATE  20,120,20,80 

,100 

640 

LAXES  5,100,0,0,2,5,5 

650 

FOR  1=2  TO  1000  STEP  2 

660 

PLOT  I*. 06, PCX) 

670 

NEXT  I 

680 

MOVE  34,0  ©  PLOT  S4,0  ©  PLOT 

690 

MOVE  S5,0  e  PLOT  S5,0  @  PLOT 

700 

MOVE  JO*. 06,0  e  PLOT  JO*. 06,0 

710 

COPY 

720 

!  FIND  M: DOT: AVE  AND  Tbar:AVE 

730 

FOR  1=1  TO  100 

740 

IF  V(I,3)<S4  THEN  770 

750 

V4  =  I 

760 

GOTO  780 

770 

NEXT  I 

780 

FOR  I=V4  TO  100 

790 

IF  V(I,3)<S5  THEN  820 

800 

V5=I 

810 

GOTO  831 

820 

NEXT  I 

830 

V6=0  e  T=0  &    T3=0 

840 

FOR  I=V4  TO  V5 

850 

V6=V6+V<I,1)  !  M:DOT  AVE 

860 

T3=T3+V<I,4)  1  T:HOT  AIR 

870 

T=T+V(I,2)  !  TjBAR  AVE 

880 

NEXT  I 

890 

V8=V5-V4  !  TIME  INT 

900 

V6=V6/V8  !  M:DOT  AVE 

910 

T3=T3/V8  1  T; HOT: AIR : AVERAGE 

920 

T=T/V8  !  T:bar  AVE 

930 

PRINT  "Pc :bar=";Pl;"  PSIA" 

940 

PRINT  '•Tbf=";S6;"  SECONDS" 

950 

PRINT  "m:AVE=";V6 

960 

PRINT  "T:BARave=";T 

970 

PRINT  "T:hot;air=";T3 

e  FRAME  e  SCALE  0,60>0 


34,P(Si) 
35,P(S2) 
e  PLOT  J0*.06,P(J0) 


91 


980  PRINT 

990  U7=W3/S6  !  FUEL  FLOW  RATE 
iOOO  Xi=U7/V6  !  fuel/air  ratio 
iOiO  PRINT  "f="iXi 

1020  IF  F1$="HTPB"  THEN  GOSUB  1130 
1030  IF  F1$="PMM"  THEN  GOSUB  1280 
1040  GOSUB  1420 
1050  GOSUB  1450 

1060  T2=(Pl*Al*C*F2/(U6+y7) )^2  I  Tcexp 
1070  N0=(T2-T)/(T1-T>  !  EFFICIENCY 
1080  PRINT  "Tc:THEORY="}Tl 
1090  PRINT  "Tc:EXP  -  '•;T2 
1100  PRINT  "T:BAR iAyE=";T 
1110  PRINT  "THERMAL  EFFICIENY  =";N0 
1120  END 

1130  !  HTPB  SUBROUTINE 
1140  F9=. 073  !  F:STOICH 
1150  D9=.0332  !  DENSITY  ♦M/in^3 
1160  X=X1/F9 

1170  T1=~933.3  +  4893*X+203*X*2-899*X^3-205.9*X''6^-1  .  1925* 
T+. 000 1014*T*2-. 0000000 184*T^3-1 . i38*T*X 
1180  T1=T1+.3864*T*X^2-.0000767*T^2*X 

1190  G=l .40162-.27723*X+. 1483*X^2-. 02448*X*3+ . 0  0523*X^6> 
-.  000  04994*T+.  0  0  00  00004255*T^2-2  .  36E-13*T''3 
1200  G=G+. 0  0OO77585*T*X-.0OOO31138*T*X^2-.0OOOO00Oi885JK 
T*2*X 

1210  R=52  .966+3  .  516#X-9  .  55*X^2+5  .  71*X''3+  .  2739*X*6-  .0008 
08*T+. 0  0  00  0  0136*T^2-5.99E-ll*T^3 

1220  R=R+.0010208*T*X-. 0001872*T*X^2+ . 0000002185*X*T^2 
1230  PRINT  "PHI=";X 
1240  PRINT  "Tc:th="jTl 
1250  PRINT  "Rc:th=^-";R 
1260  PRINT  '•GAMHAc=";G 
1270  RETURN 
1280  !  PMM  SUBROUTINE 
1290  F9=, 12  !  F:STOICH 
1300  D9=.0426  !  DENSITY  ♦M/in''3 
1310  X=X1/F9 

1320  Tl=1642-5902*X+14688*X*2-7478*X*3+4 . 05*T«X- . 18*T-4 
.68*T*X^2+1  .28*T*X*4-.  000  03259>KT*2*X^2 

1330  G=l .34-. 0265*X-.193*X"2+. 1329«X^3-. 0  00  009126*T- . 0  0 
004205*T«X+. 00007429*T*X*2-. 0000272*T*X"4 
1340  G=G+.00000000209*T*2*X"2 

1350  R=47 . 1569+X*<27 . 54+X*(2 . 9108-X*<  41 . 242-X*( 19 . 545-1 
i.567*X^3)) )) 

1360  R=R+.0  00  0  01*T*(12403-.0032029*T''2) 
1370  PRINT  "PHI=";X 
1380  PRINT  "Tc:th="jTl 
1390  PRINT  "Rc:th=";R 
1400  PRINT  "GAhMAc=";G 
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1410  RETU«N 

1420  !  F(GAhhA:C) 

1430  F2=SQR(G*G1/R*(2/(G+1) ) " ( ( G+1 ) / ( G-1 ) ) ) 

1440  RETURN 

1450  !  D: FINAL: AVERAGE 

1460  D8=SQR(4*U3/<PI*D9*L2)+D0^2)  !  D:FINAL  A^E 

1470  R8=(D8-D0)/(2*S6)  !  AVE  BURN  RATE 

1480  D8=2.54*D8  !  IN  CM. 

1490  R8=2.54«R8  !  IN  CM/SEC 

1500  PRINT 

1510  PRINT  "d:f inal!AVE=";Da 

1520  PRINT  "AyE:burn  RATE=";R8 

1530  RETURN 


93 


10  !  HELP 

20  !  BERT  HANSEN 

30  ! 

40  !  THIS  PROGRAM  INITIATES  THE  LIGHT  DIFFRACTION  EXPER 

IMENT. 

50  !  IT  ALSO  PROVIDES  HELP  WITH  THE  DOCUMENTATION. 

60  ! 

70  OPTION  BASE  i 

80  CLEAR 

90  COM  Qi$[i2] >Q2*[i2],Ti*C20] ,Di$Ct2] ,D2$[i2] 

100  DISP  "ENTER  DATE : TIME< MMDDHHMMSS ) " 

110  INPUT  Tl$  . 

120  OUTPUT  709  j"TD"iTU        .>  +  ^t,^^':v  .v.'y.f 

130  OUTPUT  709  } "TD"  ^ 

140  Q$="DATA"  !  TYPICAL  NAME 

150  CLEAR 

160  F=0  !  FLAG  RESET 

170  CRT  IS  1 

180  DISP  "ENSURE  THAT  THE  FOLLOWING  EQUIPMENT  IS  TURNED 

ON." 
190  DISP  "    l.ALL  HP-IB  EQUIP." 
200  DISP  "    2.VISIC0RDER" 
210  DISP  "    3.RETIC0N  ARRAY." 
220  DISP  "    4. LASER" 

230  DISP  "    5. BATTERY  AND  IGNITION  CONTROL  CIRCUITS." 
240  DISP  "    6. LOW  PASS  FILTER" 

250  DISP  "The  follouinq  prograMS  are  asssociated  with  t 
he  light  diffraction  experinent," 
260  DISP 

270  DISP  "    PRESS  [CONT]  WHEN  READY." 
280  IF  F=0  THEN  PAUSE 

290  DISP  "  l.ACQDTA-  acquires  data  fron  both  diode  arra 
ys  and  stores  it  on  the  disc." 
300  DISP 

310  DISP  "  2.CMBINE-coMbines  raw  data  arrays  into  TWO  1 
024x2  arrays  and  datafiles." 
320  DISP 

330  DISP  "  3.PLTDTA-plots  the  raw  data  frow  the  datafil 
e  as  VOLTAGE  VS.  DIODE  NUMBER." 
340  DISP 

350  DISP  "    PRESS  [CONT]  WHEN  READY." 
360  IF  F=0  THEN  PAUSE 

370  DISP  "  4 .REDUCE-perf orrts  the  reduction  «ath  and  plo 
ts  NORMALIZED  INTENSITY  VS.  THETA . " 
380  DISP 

390  DISP  "  5.REDUC2-plots  the  ideal  curve  fro«  disc.  Ca 
libration  data  can  also  be  plotted." 
400  DISP 

410  DISP  "    PRESS  [CONT]  WHEN  READY." 
420  IF  F=0  THEN  PAUSE 
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430  DISP  "  6.RED.UC3 -CoMputes  D32  fron  the  user  sp»icifie 

d  datafileand  calcs.  theta  bar." 

440  DISP 

450  DISP  "  7.CALIB-  If  a  run  is  to  be  a  calibration  cur 

ue  this" 

460  DISP  "  prograM  perforn^s  a  curve  fit  to  get  a  SMOOth 

curve  for  coMparison." 
470  DISP 

480  DISP  "  8.REDUC4-Plots  INTENSITY  US  THETA  BAR" 
490  DISP  "using  D32  enterred  by  user  in  fornat  like  RED 
UC2" 

500  IF  F=0  THEN  PAUSE 
SiO  IF  F=2  THEN  GOTO  580 
520  DISP  "  DO  YOU  WANT  A  HARD  COPY  OF  THE  INSTRUCTIONS? 

(Y/N)" 
530  INPUT  A* 

540  IF  A$="N"  THEN  GOTO  580 
550  F=2  !  FLAG  SET 

560  CRT  IS  2   ^- 

570  GOTO  180  !  RETURNS  TO  PRINT 

580  CRT  IS  i 

590  CLEAR 

600  CHAIN  "ACQDTA" 

610  END 
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iOO 
liO 
120 
130 
140 
150 
160 
170 
180 
190 
200 
,13 
210 
220 
230 


!  ACQDTA 

!  BERT  HANSEN 

OPTION  BASE  .1 

OUTPUT  709  "TO500"  !  TIME 

COM  Ql$[121 ,Q2$[12] ,T1$[20] , 

SHORT  A(1032) ,B(1032) 

OUTPUT  709  ; "TD" 

ENTER  709  ;  Tl$ 

PRINT  "DATE:TIME  IS 


SIGNAL 
Dl$t 


TO 
12] 


^ISICORDER 
,D2$[ 12] 


OUTPUT 
OUTPUT 
1,0,13 


709 
723 


"AR .DO4.0 


o 


Tl$ 
3" 


2,0,13 


'CC,2,3,12,13T,WF,3. 1 , 0 ,3 .2, 0 ,3 . 3 , OT, UF 
3,0T" 


CLEAR 

CLEAR  722 

DISP  "ENTER  THE  FILENAME  OF 
EATED. (RAW1,RAU2)  (MAX  12  CHAR. 
240  DISP  "DATA  SHOULD  BE  STORED 
(;D70i)  AFTER  FILENAME." 
2S0  DISP  '•  BE  SURE  THE 

DISC** 

INPUT  D1$,D2* 

CLEAR 

M5'i024 

CREATE  Di*, 1024, 16 


THE  DATA  FILES  TO  BE  CR 


) 


ON   DRIVE  *i  BY  TYPING 


FILES  NAMED  DO  NOT  ALREADY  EXIST 


ON 
260 
270 
280 
290 

300  CREATE  D2*, 1024, 16 
310  PRINT  " 
AS  FOLLOWS.  " 

320  PRINT  "NO  PARTICLES-" , Dl$ 
330  PRINT  "PARTICLES-", D2* 
340  DISP  "ENTER  THE  THRESHOLD 
ES  (PSD" 
3S0  INPUT  S8  !  DESIRED  PRESSURE 


,i«/>-^ ---  ■^^'" 


ct^ 


9 


DATA  WILL  BE  STORED  ON  DISC  WITH  FILE  NAME 


PRESSURE  TO  TRIGGER  DEVIC 


360  DISP  "ENTER  TIME 
ESSURES  (SECS. )" 
370  INPUT  T8 

T8=(T8-.25)*1000 

yO=S8/i5i.5 


DELAY  IN  SECONDS  FROM  THRESHOLD  PR 


380 

390 

400 

410 

420 

430 

440 

450 

460 

LINE 

470  LOCAL  LOCKOUT  7 

480  R7=225 

490  OUTPUT  722  "HSM002SW2SO1L1S0F4R iT3QXl " 

500  GOSUB  1260  !  ENTER  READING 

510  Ra=v 


I 

CLEAR 

DISP  "  LASER 

D3$=Di$ 

GOSUB  840  ! 

GOSUB  1180  ! 

OUTPUT  709 


ON  , IGNITION  OFF" 

CALL  MULTIPROGRAMMER 
!  STORE  DATA 
"ACO"  !  CONNECT  3456  TO 


IGNITION  FIRING 
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520  IF  R8<R7  THEN  GOTO  S70 

530  !  CONTINUITY  CHECK 

540  DISP  "CONTINUITY  CHECKS  BAD  RECHECK  BEFORE  PROCEEDI 

NG,  THEN  PRESS[CONTJ  " 

550  PAUSE 

560  GOTO  490 

570  CLEAR 

S80  DISP  "CONTINUITY  IS  OK" 

590  OUTPUT  709  ;"DC4,2"  !  NITROGEN  ON 

600  DISP 

610  OUTPUT  709  j "AC2"  !  CONNECT  D^M  TO  FIRE  SWITCH 

620  DISP  '•  ENSURE  VISICORDER  IS  SET  UP  TO  RUN  ON  PROPErt 

SCALE  AND  LAMP  IS  ON" 
b30  OUTPUT  722  "HSM002SW2Z0SOiLiFL0 . 0 iSTISOFlR iT3Q" 
640  DISP  "STANDING  BY  FOR  IGNITION" 
650  BEEP 
660  OUTPUT  722  "Xi"  ©  GOSUB  1260  !  ENTER  DATA 

670  R9=ABS(V)  6  IF  R9<iO  THEN  GOTO  660 

680  O0=TIME  8  OUTPUT  709  "DC4>0"  @    ON  TIMER*  1,5000  GOS 

UB  1250  — 

690  OUTPUT  709  "ACl"  @    CLEAR 

700  OUTPUT  722  ;"X1"  S  GOSUB  1260  !  ENTER  DATA 

710  R9-ABS(V)  e  IF  R9<<v'0  THEN  GOTO  700 
720  WAIT  T8  9    Ol=TIME 
730  GOSUB  840  !  MULTIPROGRAMMER 
740  WAIT  3000 
750  OUTPUT  709  ;"DO4,0" 
760  LOCAL  7 
770  D3$=D2* 

780  GOSUB  1180  !  STORE  DATA 

790  PRINT  "TIME  FIRE  TO  MULTI  CALL=";Oi-00 
800  PRINT  "TIME  FIRE  TO  MEM  INTR=";O2-O0 

810  PRINT  "DATA  STORED  ON  DISK  WITH  FILENAMES  i  " ',V>i% ',  * ) 
AND  (";D2*j")" 
820  CHAIN  "CMBINE" 
830  END 

840  !  MULTIPROGRAMMER  ROUTINE 

850  OUTPUT  723  ;  "SF , 2 > 3 , 1 ,  . 001 , i2T ,SF , 12 ,3 , i ,  . 00 1 , i2T ' 
860  OUTPUT  723  ; "WF ,3 . 1 , OT , WF ,3 . 2 , OT, WF ,3 . 3 , OT,UF , 13 . 1 , 
0T,WF,13.2,0T,WF,13.3,0T" 

870  OUTPUT  723  ; "WF ,3, 1031T , WF , 13 , 1031T ,MI ,2 , 1 032T" 
880  OUTPUT  709  j"DC4,3"  !  MARKS  THE  ^ISICORDER  TRACE 

890  OUTPUT  723  j  "  AC  ,3T , WF , 4 . 2 , IT , OP  ,  4 , 16 . 5T " 

900  K=SP0LL(723)  !  WAIT  FOR  MEM  INTERUPT 

910  IF  K064  THEN  GOTO  900 

920  02=TIME 

930  OUTPUT  709  i"D04,3" 

940  DISP  "  DATA  TAKEN  AND  STORED' 
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950  SEND  7  ;  UNL  MLA  TALK  23  SCG  t;^ 

960  ENTER  7  ;  Ki 

970  DISP  "MEMORY  BOARD  HAS  DATA"        <0 

980  DISP  "  SLOT*="jKi  ~^     \ 

990  OUTPUT  723  ;  "DC  ,3T , WF , 4 . 2 , OT"  ' 

1000  ENABLE  INTR  7;8 

iOiO  SEND  7  ;  UNL  MLA  TALK  23  SCG  5   _  ?'   _//, 

1020  FOR  1  =  1  TO  1025  STEP  8       _^4_/ -  !  "^r.^, - 

1030  ENTER  723  USING  "*  ,K  ,K  ,K  ,K  ,K  ,K  ,K ,K "  ;  A(I)>A<I  +  i), 

A(I+2)  ,A(I  +  3)  ,A(I  +  4)  ,A<I+5)  ,A<I+6)  ,A(I-»-7) 

1040  NEXT  I 

1050  ENTER  1 

1060  DISP  "  ARRAY  ONE  DATA  ENTERRED" 

1070  OUTPUT  723  ; "MI , 12 > 1032T" 

1080  SEND  7  ;  UNL  MLA  TALK  23  SCG  5 

1090  FOR  1=1  TO  1025  STEP  8 

1100  ENTER  723  USING  "* , K , K  J<  ,K ,K ,K ,K ,K "  )  B(I),B<I  +  1), 

B(I+2) >B(I+3) ,B<I+4) ,B(I+5) ,B(I+6),B<I+7) 

1110  NEXT  I 

1120  ENTER  7 

1130  DISP  "  ARRAY  TWO  DATA  ENTERRED" 

1140  HAT  A=-A 

1150  MAT  B=-B 

1160  OUTPUT  723  j "CC ,2 , 3 , 12 > 13T " 

1170  RETURN 

1180  !  THIS  ROUTINE  STORES  THE  RAW  DATA  ON  THE  DISK 

1190  ASSIGN*  1  TO  D3$  j  "•  ^_   -  ao' 

1200  FOR  1=7  TO  1027  STEP  4  -^ ~ 

1210  PRINT*  1  ',    A(I)  ,A(I  +  i)  ,A(I+2)  ,A<I+3)  ,B<I)  ,B(I  +  1),B  \ 

(I+2>,B(I+3)  / 

1220  NEXT  I 

1230  ASSIGN*  1  TO  « 

1240  RETURN 

1250  OFF  TIMER*  1  &    OUTPUT  709  ;"DO4,0,2"  @  RETURN 

1260  !  SRQ  SERVICE 

1270  STATUS  7,1  ;  AO 

1280  ENTER  722  ;  V 

1290  ENABLE  INTR  7j8  @  RETURN 
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10  !  CMBINE 

20  !  COMBINES  RAW  DATA  FILES  INTO  A  SINGLE  DATA  FILE 

30  OPTION  BASE  i 

40  COM  Qi*C12] ,Q2$[12] ,Ti*C20] ,Di$[i2] >D2*[i21 

50  D9$=": 0701" 

60  CLEAR 

70  ! 

80  SHORT  Y(1024,4) 

90  D3$=Di$ 

100  FOR  M=i  TO  2 

110  ASSIGN*  2  TO  D3t>  

120  FOR  1=1  TO  1021  STEP  4 

130  READ*  2  ;  Y(I,M),Y(I+1,M) ,Y(I+2,M)>Y(I+3,M) ,Y(I,M+2 

)  ,Y(I-H,M+2)  ,Y(I+2,M  +  2)  ,Y(I+3,M+2) 


ASSIGN*  2  TO  ^4^-     ^^  ..  „,,_.,X^^-f  .^T^^T^^^^ 


150 

160  D3$=D2$ 
170  NEXT  M 
180  Q3*=D1$ 

190  FOR  M=l  TO  3  STEP  2 

200  ASSIGN*  i  TO  Q3$  i^l    r' 

210  FOR  N=l  TO  1017  STEP  8 

220  A=«Y(N,M)  e  B=Y(N-f-l,M)  9  C=Y<N+2,M)  B  D=Y(N+3,M) 
230  E=Y(N+4,M)  e  F=Y(N+5,M)  8  G=Y(N+6,M)  8  H=Y(N+7,M) 
240  R=Y<N,M+i)  9   S=Y<N+i,M+l)  9    T=Y(N+2,M+i)  S  U=Y<N+3, 
M+i)  e  V=Y(N+4,M+i) 

250  W=Y<N+S,M+1)  ©  X=Y<N+6,M+1)  &    Z=Y(N+7,M+i) 
260  PRINT*  1  ;  A  ,B,C  ,D  ,E,F,G,H,R  ,  S  ,T  ,U,^  ,  W,X  ,  i'. 
270  NEXT  N 

280  ASSIGN*  1  TO  *  -^^(^ /Utxt  Ao^-^.^^  > 
290  Q3*=D2$ 
300  NEXT  M 

310  DISP  '•  DATA  STORED  ON  DISK  ON  FILE  (";Di*;")  AND  <" 
;D2*j ")" 
320  PRINT 

330  PRINT  "DATA  FILES  HAVE  BEEN  REARRANGED  A  FOLLOWS:" 
340  PRINT  "PHOTO  ARRAY  *i  IS  ";D1$ 
350  PRINT  "PHOTO  ARRAY  #2  IS  " ; D2$ 
360  CHAIN  "PLTDTA" 
370  END 

--."f)   -  ., 
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iO  •  PLTDTA   RE^JISION  i 

20  ! 

30  !  PLOTTING  RAW  DATA  FROM  ARRAY 

40  ! 

SO  !  BERT  HANSEN 

60  CLEAR 

70  OPTION  BASE  1 

80  COM  Qi$Ci2] ,Q2S[i2],Ti*[201,Di$[i23,D2$[i2] 

90  OUTPUT  709  ; "TD" 

100  ENTER  709  ;  Ti$ 

110  DISP  "  ENSURE  PLOTTER  IS  READY  FOR  PLOTTING  WITH  CL 

EAN  PAPER  AND  DESIRED  PEN." 

120  DISP  "PRESS  [CONT]  WHEN  READY!" 

130  BEEP 

140  PAUSE 

150  DISP  "  ENTER  THE  FILE  NAME  OF  THE  DATA  TO  BE  PLOTTE 

D.  (";D1*;  "  OR  ";D2»; "  )" 

160  DISP  "DRIVE  MUST  BE  SPECIFIED  AS  (NAME:D701)" 

170  INPUT  Q3$ 

180  DISP  "ENTER  DESIRED  PLOTTING  INTERVAL  OF  DIODE  ARRA 

Y  (2,4,6)" 

190  INPUT  J 

200  CLEAR 

210  DISP  "READING  DATA  FROM  DATA      FILE  ON  DISC  DRIVE 

M 

220  PLOTTER  IS  705 
230  SHORT  Y(i024,2) ,X(1024> 
240  M5=1024 

250  REDIM  X<M5) >Y(M5>2) 
260  Pli="VOLTAGE" 
270  P2$="DI0DE  NUMBER" 
280  P3$="V0LTAGE  VS.  DIODE" 
290  P4$="PLTDTA  RESULTS" 
300  ! 

310  ASSIGN*  1  TO  Q3* 
320  FOR  1=1  TO  M5-7  STEP  8 

330  READ*  1  ;  A  ,  B  ,C  ,D  ,E  ,F  ,G  ,H  ,R  ,5  ,  T  ,U  ,  V  ,  W  >  X  ,  Z 
340  Y(I,1)=A  9  Y<I-H,1)=B  ©  Y(I+2,1)=C  9    Y(I+3,i)=D  &    Y 
(I+4,1)=E  e  Y(I+5,1)=F 
350  Y<I+6,i)==G  e  Y<I+7,1)=H 

360  Y(I,2)=R  9    Y(I+1,2)=S  @  Y(I+2,2)=T  9    Y<I+3>2)=U  9   Y 
(I+4,2)=V  9    Y<I+5,2)=W 
370  Y(I-»-6,2)=X  e  Y(I+7,2)-Z 
380  NEXT  I 
390  M7=MAXAB(Y> 

400  DISP  "DATA  INPUT  COMPLETE" 
410  LIMIT  0,250,0,180 
420  FRAME 

430  LOCATE  20,120,20,80 
440  FRAME 

100 


4S0 

SCALE 

0,M5,-. 1,M7+. 1 

460 

CSIZE 

3,  .5,0 

470 

DEG 

480 

LDIR 

Q 

490 

FXD  0 

>2 

500 

LAXES 

M5/20, .02,0,0,2,5,5 

510 

MOVE 

.65*M5, .93*(M7+. 1) 

520 

LABEL 

P4$ 

530 

MOVE 

.65*M5,  .84*(M7+. 1) 

540 

LABEL 

P3$ 

550 

LABEL 

Tl$ 

560 

MOVE  1 

M5/2,-.  OV:", 

570 

LABEL 

P2$ 

580 

LDIR 

90 

590 

MOVE  - 

-< .i*M5), .5*M7 

600 

LABEL 

Pl$ 

610 

LDIR  1 

[) 

620 

M  =  l 

630 

GOSUB 

760  I  PLOT  FIRST  ARRAY 

640 

DISP 

"  FIRST  ARRAY  PLOTTED!" 

650 

M=2 

660 

GOSUB 

760  !  PLOT  SECOND  ARRA49 

2  MOVE  . 

+  .  1) 

670 

DISP 

"  SECOND  ARRAY  PLOTTED!" 

680 

DISP 

»  PLOT  COMPLETE!  * 

690 

DISP 

"DO  YOU  WISH  TO  PLOT  THE 

OTHER  PH 

E? 

<Y/N)" 

700 

INPUT 

N$ 

710 

IF  N*= 

="Y"  THEN  GOTO  110 

720 

CLEAR 

730 

DISP 

"LOAD  AND  RUN  < REDUCE)  TO 

PROCEED 

DUCTIONS. 

n 

750 

END 

760 

!  PLOT  DATA 

770 

X(l)  = 

1 

780 

MOVE 

X(1),Y(1,1> 

790 

FOR  I 

=1  TO  MS  STEP  J 

800 

X(I)  = 

I 

310 

PLOT 

X<I),Y(I,M> 

820 

NEXT 

I 

830 

PENUP 

840 

RETURN 

.65*M5, .74*(M7 


ARRAY  FIL 


WITH  DATA  RE 


101 


10  I  REDUCE 

?.0  CLEAR  7 

30  DISP  "REDUCE  LOADED" 

40  ! 

SO  !  BERT  HANSEN 

bO  !  THIS  PROGRAM  REDUCES  THE  ARRAY  DATA 

70  OPTION  BASE  i 

80  DIM  A(2,2) ,B(2) ,Yi(2) 

90  COM  Q3$[i2] ,Q2$[i2] >Ti$[20] ,Dl*[i21,D2$[i2] ,« 

100  OUTPUT  709  ; "TD" 

110  ENTER  709  ;  Tl* 

120  CLEAR 

130  DISP  "DATA  REDUCTION  BEGINS" 

140  DISP  "  ENTER  THE  INDEX  OF  REFRACTION^? < AIR  M=i.O,  W 

ATER  M=i.3S)" 

ISO  INPUT  M 

160  DISP  "INSTALL  PAPER  AND  PEN  ON  PLOTTER  TO  PLOT  I() 

VS  '• 
170  DISP  -PRESSCCONT]  WHEN  DONE" 
180  PAUSE 

190  DISP  "ENTER  THE  FILE  NAME  OF  DATA  TO  BE  REDUCED." 
200  DISP  "INCLUDE  THE  DISC  LOCATION  ON  FILENAME  (  NAME: 
D701)" 

210  INPUT  A^:| 
220  D9$="!D701" 

230  L=:.  000000488  !  WAVELENGTH  IN  METERS 
240  SHORT  Ti(1024) ,11(1024) ,12(1024) 
2S0  M5=1024  !  CURRENT  NUMBER  OF   ACTUAL  READINGS 
260  D8=1.5365  !  MMETERS 
270  F8=:S26  !  MMETERS 
280  Dl=.025  !  MMETERS 
28S  GCLEAR 
290  PLOTTER  IS  70S 
300  LIMIT  0,250,0,130 
310  LOCATE  20,120,20,80 
320  FRAME 

330  SCALE  0,  .04, -.1,1 
340  CSIZE  3, .5,0 
350  DEG 
360  LDIR  (1 
370  FXD  3,1 

380  LAKES  .002,-1,0,0,2,2,5 

390  DISP  "DO  YOU  WANT  TO  PLOT  A  COMPARISON  CURVE  (  UNIC 
UR)?  (Y/N)" 
400  INPUT  N* 

410  IF  N*="Y"  THEN  GOSUB  1140 
420  ASSIGN*  1  TO  A$ 
430  FOR  1=1  TO  M5-7  STEP  8 

440  READ*  1  ;  A ,B , C ,D ,E,F ,G ,H ,R , S , T ,U , V , W , X , Z 
450  I2(I)=A  e  I2(I  +  1)=B  @  I2(I-»-2)^C  e  I2(I-»-3)=D  9    12(1  + 
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4)=E  e  12<l^^j)=F    e  I2(I+6)=G  &  I2<I+7)=H 

460  Ii(I)=R  8  Ii(H-i)=S  e  Ii(I  +  2)-T  e  Ii(I+3)=U  @  Ii(I  + 

4)=v  e  ii(i+5)=u  e  ii(i+6)=^x  e  ii(i+7)^z 

470  NEXT  I  '     ,  ,  ,  V 

480  MAT  Ii  =  Ii-i;^i    /^;!;^-/r/<-      -^  ^/^/?/ '^/.Sj 

490  Z9=AMAX(Ii) 

500  Z0=AMAXRQW 

5iO  PRINT  "MAXIMUM  OUTPUT  AT  DIODE#=";Z0 

520  DISP  "ENTER  THE  INCREMENT  AFTER  PEAK  <  J=:*START  ,K==I»   9 

STOP  FIT)"  p 

530  INPUT  J,K 

540  N0=K-J 

550  Xi  =  0  e  X2=0  8  X3=0  9    Yi  =  0  S  Y2-=0 

560  FOR  N=J  TO  K 

570  Ti(N)=N+6i .46 

580  Xl=Xi+Ti(N)  '  ' 

590  X2  =  X2+T1(N)^2  ^J  <-     "?   y         /=>^   </ 

600  X3=X3  +  Ti(N)*Ii(N)      J  -       /    Y    -Q-      ^^^^^ 

6io  Yi=Yi+ii(N)         \-niAA<:  V -X -^  y>-j^    ^^ 

620  Y2=Y2+Ii(N)''2  y/^/  J 

630  NEXT  N  .  ^   ' 

640  A(i>i)=NI  ^^  '  '     '^ 

650  A(i,2)=Xi 

660  A(2,i)=X.l 

670  A(2,2)=X2 

680  Yi(i)=Yi 

690  Yl(2)='X3 

700  MAT  B=SYS(A,Yi) 

7iO  PRINT  "DELTA=";B(i) 

720  PRINT  "SLQPE=";B(2)  / 

730  A3=B<1)  n}/ 

740  R2«<N0«X3-Xi*Yi)*2/<  (N0)(cX2-Xi*2)*(N0*Y2-Yi*2)  )  !  CO 

RRELATION 

750  PRINT 

760  PRINT  "CORRELATION  R''2="}R2 

770  MAT  Ii=<i/A3)*Il 

780  !  li  ARRAY  IS  THE  FINAL  PLOTTING  NORMALIZED  ARRAY 

790  !  NOU  PLOT  li  VS.  Ti 

800  FOR  I=i  TO  M5 

810  T1(I)-(I*.025+D8)/(F8*M)  !  THETA   f^/?^//  'S  A-**-/^  T/e  v/ 

820  PLOT  Ti(I),Ii<I)  ^ 

830  NEXT  I 

840  PENUP 

850  MO^IE  0,0 

860  MOVE  .02>-.3 

870  LABEL  "THETA  (RADIANS)" 

880  MOVE  -.004,-2  »  DEC  8  LDIR  90 

890  LABEL  "NORMALIZED  INTENSITY" 

900  LDIR  0  8  MOVE  .015, .8 

910  LABEL  "INTENSITY  VS.  THETA" 
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920  hOVE  , OlS, .9 

930  LABEL  "REDUCE   RESULTS" 

940  hOUE  .015,  .7 

950  LABEL  Ti* 

960  Q5$="REDUC3" 

970  DISP  "  IS  THIS  A  CALIBRATION  RUN?  (Y/N)» 

980  INPUT  N$ 

990  IF  N$="N"  THEN  GOTO  iOiO 

1000  Q5$="CALIB" 

1010  DISP  "ENTER  THE  NAME  OF  THE  DATA  FILE  TO  STORE  THE 

DATA  (INTENSITY  VS.THETA)" 
1020  DISP  "INCLUDE  THE  DISC  LOCATION  ON  FILENAME  (  NAME 
:D701)" 

1030  INPUT  Q3$  ■ — ^. 

1040  CREATE  Q3$,1024,16 

1050  PRINT  "INTENSITY  VS .  THETA  STORED  IN  FILE=:''jQ3$ 
1060  ASSIGN*  2  TO  Q3$  ,       UyUjcf^   ^^ 

1070  FOR  1  =  1  TO  MS  ^    .       _^ 

1080  PRINT*  2  ;  TKD.IKI)  '    T  ti' ^S^  ^<^'^Jr 

1090  NEXT  I  \         '^rOi?^--    ''^f 

1100  ASSIGN*  2  TO  >i 

1110  DISP  "DATA  STORED  ON  DISK  IN  FILE  —  *';Q3$ 
1120  CHAIN  Q5$ 
1130  END 

1140  !  COMPARE  WITH  UNICUR 
1150  ASSIGN*  2  TO  "UNICUR" 
1160  FOR  1=1  TO  100  STEP  2 

1170  READ*  2  ;  Tl ( I ) , II ( I ) ,T1 ( I+l ) , 11 < I+l) 
1180  NEXT  I 
1190  DISP  "ENTER  EXPECTED  VALUE  OF  Dbar  ( 23 . 98E-6 )ETC 

1200  INPUT  D 

1210  FOR  1=1  TO  99 

1220  Tl(I)=Ti(I)*L/<PI*D)  ~^' 

1230  PLOT  Tl(I)  ,11(1) 

1240  NEXT  I 

1250  PENUP 

1260  RETURN 


D^/^ 
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10  '  CALIB 
?.Q    OPTION  BASE  i 

30  COM  Q3*[i2]  ,Q2$[12] /n$t20]  ,Dl$tl21  ,D;i$[12],rt 
40  D9$=":D701" 
50  !  CALIBRATION  PROCEDURE 
60  !  BERT  HANSEN 
70  CLEAR 
80  GCLEAR 

90  DISP  "ARE  YOU  RUNNING  CALIB  SEPARATELY?  (Y/N)" 
100  INPUT  N* 
110  IF  N$="N"  THEN  160 

120  DISP  "ENTER  DATAFILE  NAME  AND  DRIVE  ADDRESS" 
130  INPUT  Q3$ 

140  DISP  "ENTER  INDEX  OF  REFRACTION  <AIR=i.O  ,WATER=i.3 
S  >  " 

150  INPUT  H 

160  DISP  "ENTER  THE  ADDRESS  OF  PLOTTER.  <CRT=1,  7225B=7 
05)" 

170  INPUT  P 
180  PLOTTER  IS  P 

190  IF  P---705  THEN  LIMIT  0,250,0,180 
200  LOCATE  20,120,20,80 
210  FRAME 

220  SCALE  0,10,-. 1,1 
230  LAXES  .5, .2, 0,0, 2,5, S 
240  ! 

250  DIM  A(ii,il) ,B(li),W(20) ,Yi(ll) ,Y2(il) 
260  SHORT  X(1024) ,Y(1024) 
270  MAT  A=ZER 
280  MAT  B^ZER 
290  MAT  Y1=ZER 
300  MAT  W=ZER 
310  M5=i024 

320  DISP  "  INDEX  OF  REFRACTION  IS  =";M 
330  DISP  "USING  <  ">Q3»;"  >" 

340  DISP  "A  POLYNOMIAL  REGRESSION  ROUTINE  OF  (J)TH  ORDE 
R  WILL  BE  PERFORMED.  " 

350  DISP  "ENTER  ORDER  OF  POLYNOMIAL<  REASONABLE).- 
360  INPUT  Ji 
370  ASSIGN*  1  TO  Q3* 
380  FOR  1=1  TO  M5  STEP  2 
39a  READ*  1  J  X<I) ,Y(I) ,X<I+i) ,Y(I+1) 
400  NEXT  I 

410  DISP  "ENTER  Dbar   TO  USE  IN  DATA  REDUCTION  (e.g.  23 
.98E-6) " 
420  INPUT  D3 
430  L=. 000000488 
440  F=526  !  FOCAL  LENGTH 
450  K6=D3*PI/L 
460  FOR  1=1  TO  MS 
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470  X(I)=X(I)*K6 

480  PLOT  Xvl) ,Y(1) 

490  NEXT  I 

500  !  I  VS.THETA  BAR 

SIO  PENUP 

S20  DISP  "ENTER  THE  MINIHUH  THETA  BAR  OF  INTEREST." 

530  INPUT  J 

540  J=M*J/(236*D3)-6t .46 

550  K=i .5365 

560  DISP  "ENTER  THE  MAXIMUM  THETA  BAR  OF  INTEREST  .(ba 

r=3.5)" 

570  INPUT  Mi 

580  Mi=(Ml*L*F*M/(PI*D3)-K)/. 0  2533 

590  PRINT  "MIN  DIODE  =-;J 

600  PRINT  "MAX  DIODE  ="jMi 

6iO  J2=Ji+i 

620  FOR  I=i  TO  2*Ji 

630  ! 

640  S0=0 

650  FOR  M=J  TO  Mi 

660  ! 

670  SO=X(M)^I+SO 

680  NEXT  M  ~^ 

690  W(I)=SO 

700  NEXT  I 

7iO  ! 

720  FOR  I=i  TO  J2 

730  S0=0 

740  FOR  N=J  TO  Mi 

750  SO=SO+Y(N)*X(N)^<I-i) 

760  NEXT  N 

770  Yi(I)=SO 

780  NEXT  I 

790  NO=Mi-J  !  *  OF  DATA  PAIRS 

800  ! 

8iO  REDIM  A<J2,J2) ,B<J2) ,W<2*Ji),Yi(J2),Y2<J2) 

820  A(i,i)=NO 

830  FOR  1=2  TO  J2  !  FIRST  ROW 

840  A<i,I)=W(I-i) 

850  NEXT  I 

860  ! 

370  FOR  1=2  TO  J2  !  FIRST  COLUMN 

880  A<I,i)=W<I-i> 

890  NEXT  I 

900  I 

9iO  M=0 

920  FOR  1=2  TO  J2  •  CORE  iO*iO 

930  FOR  K=2  TO  J2 

940  A(I,K)=U(K+M) 

950  NEXT  K 

960  h=M+i 
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970  NEXT  I 

980  GOSUB  1300  !  CORRELATION 

990  ! 

iOOO  MAT  B=SYS<A,Yi) 

iOiO  PRINT 

1020  PRINT  "POLYNOMIAL  COEFFICIENTS^ : " 

1030  PRINT 

1040  MAT  PRINT  B 

1050  I 

1060  !  COMPUTE  THE  CALIBRATION  CURUE 

1070  S=.l 

1080  FOR  1=1  TO  100 

1090  X<I)=I*S 

1100  30=0 

1110  FOR  N=i  TO  JZ 

1120  SO  =  SO+B(N)*X<I)''(N-i) 

1130  NEXT  N 

1140  Y<I)=S0 

1150  PLOT  X(I) ,Y(I) 

1160  NEXT  I 

1170  DISP  "  COMPLETE  ! " 

1180  DISP  -ENTER  THE  DESIRED  FILE  NAME  TO  STORE  DATA 

1190  INPUT  Q3$ 

1200  CREATE  Q3$,100,16 

1210  ASSIGN*  2  TO  Q3$ 

1220  !  STORE  SMOOTH  CAL  CURVE  ON  DISC. 

1230  FOR  1=1  TO  100  STEP  2 

1240  PRINT*  2  ;  X ( I ) , Y ( I ) , X < I+l ) , Y ( I+l ) 

1250  NEXT  I 

1260  DISP  "DATA  STORED  ON  DISC" 

1270  ASSIGN*  2  TO  * 

1280  CHAIN  ••REDUC2" 

1290  END 

1300  !  COMPUTE  CORRELATION  FACTOR 

1310  DISP  "  CORRELATION  " 

1320  SHORT  X4(1024)  ,Y4(1024) 

1330  MAT  X4-ZER 

1340  MAT  Y4=ZER 

1350  MAT  X4-X<J:Mi) 

1360  MAT  Y4=Y(J:M1) 

1370  Ci=D0T<X4,X4) 

1380  C2=D0T(Y4,Y4) 
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1390  C3=D0T(X4,Y4) 

1400  C4=Ci^2 

1410  C5^C2*2 

1420  C6=SUM(X4) 

1430  C7=SUM(Y4) 

1440  N2=N0*C3-C6*C7 

1450  D2=SQR( (N0*C1-C6"2) *(N0»C2-C7"2  )  > 

1460  R=N2/D2  !  CORRELATION  FACTOR 

1470  R2=R^2 

1480  PRINT  "  CORRELATION  FACTOR" 

1490  PRINT  ''R^2=";R2 

1500  RETURN 
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iO  I  REDUC2 

20  ! 

30  !  BERT  HANSEN 

40  !  THIS  PROGRAM  WILL  RETREI^E   THE  UNI'^ERSAL  CURVE  DA 

TA  FROM  THE  DISC. 

SO  !  IT  WILL  THEN  PLOT  THE  CURVE  AND  THEN  RETREIVE  AND 

PLOT  THE   CALIBRATION  CURVES. 

60  ! 

70  CLEAR 

80  OPTION  BASE  1 

90  COM  Q3$[i2] ,Q2$[121 ,Ti$[20] ,Di$[ 121 ,D2$[i2] ,M 

100  D9*="!D70i" 

iiO  L=. 000000488 

120  F=526 

130  Kl=i .5365 

140  PLOTTER  IS  70^i 

150  CSIZE  3, .4,0 

160  OUTPUT  709  ; "TD" 

170  ENTER  709  j  Ti* 

180  SHORT  X(iOO),Y<iOO) 

190  SHORT  T<1024) ,1(1024) 

200  P4$="THETA  BAR " 

210  P5*=s"INTENSITY" 

220  DIM  P3*C26] 

230  P3*="INTENSITY  VS.  THETA  BAR" 

240  LIMIT  0,250,0,180 

260  LOCATE  20,120,20,80 

270  FRAME 

280  SCALE  0,10,0,3 

290  LINETYPE  1 

300  E=i 

310  E3=3 

320  E1$="A,.3D-  e  E2$=".3D,A- 

330  Y=i  ©  Y2=iO 

340  GOSUB  1230  !  LABEL 

350  GOSUB  1300  !  LOG  SCALE 

360  E1$="A,.2D"  e  E2$=".2D,A" 

370  E=iO*E 

380  Y=20  9    Y2=100 

390  GOSUB  1300 

400  E1*="A,D.D"  ©  E2*="D.D,A" 

410  Y=200  ©  Y2=1000 

420  E=10*E 

430  GOSUB  130  1 

440  DEG 

450  LDIR  -91 

460  FOR  X=0  TO  10 

470  !  YAXIS  X 

480  MOVE  X,-.l 

490  LABEL  USING  "DD"  ;  X 
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50  0  NEXT  X 

5iO  LDIR  (I 

520  DISP  "SCALE  COMPLETE" 

530  SCALE  0,10,0>i 

540  CSIZE  4, .5,0 

550  FXD  1,2 

560  MOVE  5,-. 15 

570  LABEL  P4$ 

580  DEG 

590  LDIR  90 

600  MOUE  -1, .5 

610  LABEL  P5* 

620  LDIR  0 

630  MOVE  7, .8S 

640  !  LABEL  P3$ 

650  MOVE  7, .95 

660  !  LABEL  "REDUC2  RESULTS" 

670  MOVE  7, .75 

680  !  LABEL  Tl* 

690  Q5$="UNICUR" 

700  DISP  "  LABELING  COMPLETE" 

710  DISP  "CHANGE  PENS  AND  LCGNTr' 

720  PAUSE 

730  SCALE  0,10,-3,0 

740  M5=«i024 

750  ASSIGN*  1  TO  Q5$ 

760  FOR  1=1  TO  99 

770  READ*  1  ',    X(I),Y(I) 

780  PLOT  X(I) ,LGT(Y(I)) 

790  NEXT  I 

800  PENUP 

810  ASSIGN*  1  TO  * 

820  DISP  "  DO  YOU  WANT  TO  PLOT  A  CALIBRATION  CURVE?  <Y/ 

N)" 

830  INPUT  N$ 

840  IF  N*="N"  THEN  GOSUB  930 

850  DISP  "  ENTER  THE  FILE  NAME  OF  THE  CALIBRATION  CURVE!: 

TO  BE  PLOTTED.  (CALi,ETC)" 
B60  DISP  ••  INCLUDE  THE  DISC  DRIVE  WHERE  CAL  FILE  IS  STO 
RED — E.G.  -NAME:D701" 
870  INPUT  Q3$ 
880  !  DISP  "ENTER  THE  NUMBER  OF  DATA  POINTS  IN  THE  FILE 

. (100, 1024, ETC)" 
890  !  INPUT  M2  !  CAL  FILES  HAVE  100  ELEMENTS 


!  CALC  THETA  BAR  AND  PLOT  ! 


900 

M2=100 

910 

GOSUB  980 

920 

GOTO  820 

930 

LINETYPE  1 

940 

CLEAR 
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950  DISP  "PLOTS  COMPLEfE  !" 

960  DISP  "LOAD  REDUC3  TO  COHPUTE  D32  FROM  A  DATA  SET." 

970  END 

980  !  READ  AND  PLOT  CAL 

990  DISP  "ENTER  THE  INDEX  OF  REFRACTION  USED  FOR  THE  CA 

L  RUN.  (AIR  =  i.O,  WATER-=i  .35)" 

iOOO  INPUT  H 

iOiO  DISP  "ENTER  THE  Dbar  USED  TO  COMPUTE  THETA  BAR." 

1020  INPUT  D 

1030  DISP  "ENTER  THE  DIODE  NUMBERS  FOR  THE  REGION  OF  IN 

TEREST.  (P1,P2)  TYPICALLY  (50,400)'* 

1040  INPUT  P1,P2 

1050  H9=D*PI«.025/(L*F*M> 

1060  Pl=H9*(6t .46+Pl) 

1070  P2=H9*(61  .46-»P2) 

1080  PRINT  "LOWER  LIMIT  OF  THETA  BAR=";P1 

1090  PRINT 

1100  PRINT  "UPPER  LIMIT  OF  THETA  BAR=";P2 

1110  ASSIGN*  1  TO  Q3$ 

1120  FOR  1=1  TO  M2 

1130  READ*  1  ;  T(I) ,I< I> 

1140  NEXT  I 

1150  MOVE  T(1),I( 1) 

1160  FOR  1=1  TO  100 

1170  IF  TdXPi  THEN  1200 

1180  IF  T(I)>P2  THEN  1200 

1190  IF  I<I)>0  THEN  PLOT  T ( I ) ,LGT ( I ( I > ) 

1200  NEXT  I 

1210  PENUP 

1220  RETURN 

1230  !  LABEL 

1240  Y3=Y/1000 

1250  MOVE  10.i,LGT<Y)  8  LORG  2 

1260  LABEL  USING  El*  ;  "-'•,Y3 

1270  MOVE  -.i,LGT(Y)  &    LORG  8 

1280  LABEL  USING  E2$  ;  Y3,"-* 

1290  RETURN 

1300  FOR  Y=Y  TO  Y2  STEP  E 

1310  E2=LGT<Y) 

1320  !  XAXIS  LGT(Y) 

1330  IF  RMD<Y,10*E)*0  THEN  1350 

1340  GOSUB  1230  !  LABEL 

1350  NEXT  Y 

1360  RETURN 
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10  !  REDUC3 
40  ! 

50  !  THIS  ROUTINE  COhPUTE  THE    ^^ALUE  OF  D32  FOR  THE  DA 
TA  IN  THE  FILE  "Q$- 
60  !  BERT  HANSEN  13  JULY  82 
70  ! 

80  OPTION  BASE  i 
90  M5=i024 

iOO  COM  Q3$[i2] ,Q2$[12] ,Ti$[20] ,Di$[i2] ,D2$tl2] ,M 
110  CLEAR 
120  D9$="  :I)701'• 
130  SHORT  R(20) ,Xi(100)  ,Y1(100) 

140  SHORT  X(1024) ,Y(1024) ,D3<20) /r3(20) ,T(1024) 
150  DISP  "ENTER  THE  FILENAME  OF  THE  DATA  TO  BE  REDUCED 
(DATA! :D70i,ETC)" 
160  INPUT  Q3* 
170  ASSIGN*  2  TO  Q3i$ 
180  FOR  1=1  TO  M5-3  STEP  4 

190  READ*  2  ',    X<I)  ,Y(I),X(I  +  i)  >Y<I  +  1)  ,X(I+2>  ,Y(I  +  2)  ,X(I 
+3),Y(I+3) 
200  NEXT  I 

210  !  FOR  1=1  TO  1024 
220  !  X(I)=X(I)/198.5<i 
230  !  NEXT  I 

240  !  THETA  AND  INTENSITY  ENTERRED 
250  L^". 000000488 
260  Ki=L/PI 

270  DISP  "ENTER  THE  NAME  OF  THE  DATA  FILE  OF  THE  REFERE 
NCE  ARRAY. ( UNICUR ,CALi ,CAL2, CAL3) " 
280  INPUT  Qi* 
290  M4=:iOO 

300  ASSIGN*  1  TO  Ql$ 
310  FOR  1=1  TO  M4-i  STEP  2 

320  READ*  1  ;  XI ( I ) , Yi < I ) ,Xi ( I+i ) , Yi ( I+l ) 
330  NEXT  I 

340  DATA  .5, .3, .2, . 1, . 09, .08, .07, .06, .05, .04, . 03, .02, , 0 
1,  .008, .005 
350  FOR  N=i  TO  15 
360  READ  R(N)  !  REF  INTENSITIEii 
370  NEXT  N 

380  DISP  "ENTER  THE  NUMBER  OF  DIODES  TO  BE  OMMITTED  CLO 
SE  TO  THE  CENTERLINE.  " 
390  INPUT  J 
400  !  D32  CALCULATION 
410  Jl=2 

420  FOR  D=i  TO  15 
430  Ri=R<D)  !  REF  INTENSITY 
440  01=AMIN(Y)  !  SMALLEST  INTENSITY 
450  IF  RKOi  THEN  GOTO  1250 
460  FOR  I=Ji  TO  100 
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470  IF  Yi(I)>Rl  THEN  GOTO  610 

480  IF  Yl(I)>=Yi(I~i)  THEN  GOTO  600 

490  !  FIND  THETA  BAR 

500  H2=Yi(I-t) -Yi(I) 

510  K2=Xi(I)-Xi(I-i) 

520  H3=Yi(I-i)-Ri 

530  K3-K2*(H3/H2) 

540  T3(D)=«Xi(I-i)+K3  •  THETA  BAR 

550  PRINT  •'I  =  *';I;Yi(I)  jRi;D 

560  PRINT  "I  =  '*;I;T3<D)  iXi(I-i)+K3 

570  PRINT 

580  J1=I 

590  GOTO  620 

600  PRINT  "  ZERO  OR  POS  SLOPE" 

610  NEXT  I 

620  !  FIND   THE  INTENSITY  CORRESSPONDING  TO  DESIRED  INT 

ENSITY 

630  h9=i000 

640  L9=i 

650  FOR  I=J  TO  H9 

660  IF  Yd)  >0  THEN  691 

670  PRINT  "INTENSITY  VECTOR  GOING  NEGATIVE  AT  I=";I 

680  PRINT  "INCREASING  J=J+5"  @  J=J+5  8  GOTO  630 

690  Z=0 

700  FOR  M=I-5  TO  1+5 

710  Z*Z+Y<M) 

720  NEXT  H 

730  Z:=.i*2 

740  IF  Z>Ri  THEN  GOTO  790 

750  !   Z  IS  LESS  THAN  Ri 

760  J=I  !  STARTING  POS  FOR  NEXT  SEARCH 

770  DISP  "  I  FOUND  IT" 

780  GOTO  810 

790  L9'=Z  !  LAST  AVERAGE!: 

800  NEXT  I 

810  !  FIND  THETr^ 

820  K7=X(J)-X(J-i) 

830  H7=L9-Z 

840  H6=L9-Ri 

850  K6=K7*(H6/H7> 

860  T4=X<J-1)+K6  !  THETA 

870  !  COHPUTE  D32 

880  D3(D)=Ki*T3(D)/T4  !  D32 

890  PRINT 

900  PRINT  "FOR  INTENSITY^ " ; R < D) 

910  PRINT  "THETA=  ";T4;"  RADIANS" 

920  PRINT  "THETA  BAR=  "iT3(D) 

930  PRINT  "D32=   " ; D3(D>*i0OOOO0 ; "   MICRONS" 

940  NEXT  D 

950  D5=0 
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960  FOR  Ii=i  TO  D 

970  D5=DS+D3(Ii) 

980  NEXT  li 

990  D5=D5/D*iOOUOO() 

iOOO  PRINT  "032  A^^ERACE  =  "  ;  DS ;  "  MICRONS  - 

iOiO  PRINT 

1020  DISP  "DO  YOU  WANT  TO  COMPARE  DATA  TO  ANOTHER  REFER 

ENCE  FILE?  (Y/N)" 

1030  INPUT  N$ 

1040  IF  N$="N"  THEN  1080 

1050  DISP  "ENTER  THE  FILENAME  OF  THE  REF .  DATAFILE  (CAL 

1A,ETC)" 

1060  INPUT  Ql$ 

1070  GOTO  290 

1080  DISP  "DO  YOU  WANT  TO  REDUCE  ANOTHER  DATAFILE?  (Y/N 

)" 

1090  INPUT  H% 

1100  IF  N$="N"  THEN  1180 

1110  DISP  "ENTER  THE  FILENAME  OF  THE   DATAFILE  TO  BE  RE 

DUCED.  (DATA!, ETC)" 

1120  INPUT  Q* 

1130  ASSIGN*  1  TO  Q$ 

1140  FOR  I»l  TO  1021  STEP  4 

1150  READ*  1  }  X(I),Y(I)  ,X(I  +  1)  ,Y(IH)  ,X<I+2),Y<I+2),X( 

1+3), Y< 1+3) 

1160  NEXT  I 

1170  GOTO  400  !  REPEAT  D32  CALCULATIONS 

1180  DISP  "REDUCTIONS  COMPLETED" 

1190  DISP  "WANT  TO  DO  ANOTHER  RUN?" 

1200  INPUT  N$ 

1210  IF  N*="N"  THEN  1230 

1220  GOTO  10 

1230  CHAIN  "REDUC4" 

1240  END 

1250  !  REF  INTENSITY  IS  LESS  THEN  SMALLEST  INTENSITY 

1260  !  DON'T  SEARCH  JUST  COMPUTE  D32  AVERAGE  TO  THIS  PC) 

INT 

1270  D5=0 

1280  FOR  11=1  TO  D~l 

1290  D5=D5+D3<I1) 

1300  NEXT  11 

1310  D5=D5/(D-1)  !  D32  A^^ERAGE 

1320  GOTO  1000 
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iO  !  REDUC4 

20  ! 

30  !  BERT  HANSEN 

40  '  THIS  PROGRAM  WILL  RETREIVE   THE  UNIVERSAL  CURVE  DA 

TA  FROM  THE  DISC. 

50  !  IT  WILL  THEN  PLOT  THE  UNIVERSAL  CURVE. 

60  !  THE  USER  MAY  THEN  PLOT  INTENSITY  VERSUS  THETA  BAR 

WITH  D32  ENTERRED  FROM  KEYBOARD. 

70  ! 

80  CLEAR 

90  GCLEAR 

100  OPTION  BASE  t 

liO  COM  Q3$[i2],Q2$[i2],Ti*C203,Di$[i21 ,D2$[i2],M 

120  D9$=" :D70i" 

130  L=. 000000488 

140  F=526 

ISO  Ki=i.536S 

160  DISP  "ENTER  THE  PLOTTER  ADDRESS  (CRT  =  i,  7225B==705>" 

170  INPUT  P 

180  PLOTTER  IS  P 

190  IF  P=»705  THEN  LIMIT  0,250,0,180 

200  CSI2E  3, .4,0 

210  OUTPUT  709  ;"TD" 

220  ENTER  709  j  Tl$ 

230  SHORT  X(iOO) ,Y(100) 

240  SHORT  T(1024) ,1(1024) 

250  P4*="THETA  BAR " 

260  P5$="INTENSITY" 

270  DIM  P3$C26] 

280  P3$=-INTENSITY  VS.  THETA  BAR" 

290  FRAME 

300  LOCATE  20,120,20,80 

310  FRAME 

320  SCALE  0,10,0,3 

330  LINETYPE  1 

340  E=l 

350  E3=3 

360  Ei$="A,.3D"  e  E2*=".3D,A" 

370  Y=i  S  Y2=iO 

380  GOSUB  1260  !  LABEL 

390  GOSUB  1330  !  LOG  SCALIE 

400  E1*="A,.2D"  8  E2*=".2D,A" 

410  E=10*E 

420  Y=20  8  Y2=i00 

430  GOSUB  1330 

440  E1$="A,D.D"  e  E2$="D.D,A" 
450  Y=200  9  Y2=iOOO 
460  E=10*E 
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470  GOSUB  1330 

480  DEG 

490  LDIR  -91 

500  FOR  X=0  TO  10 

510  YAXIS  X 

520  HOVE  X,-. 1 

530  LABEL  USING  "DD"  ;  X 

540  NEXT  X 

550  LDIR  I 

560  DISP  "SCALE  COMPLETE" 

570  SCALE  0,10,0,1 

580  CSIZE  4, .5,0 

590  FXD  1,2 

600  MOVE  5,-. 15 

610  LABEL  P4$ 

620  DEG 

630  LDIR  91 

640  MOVE  -1,  .5 

650  LABEL  P5* 

660  LDIR  0 

670  MOVE  7,  .aSi 

680  LABEL  P3* 

690  MOVE  7, .95 

700  LABEL  "REDUC4  RESULTS" 

710  MOVE  7, .75 

720  LABEL  Ti$ 

730  Q5$='"UNICUR" 

740  DISP  "  LABELING  COMPLETE" 

750  DISP  "CHANGE  PENS  AND  IC0NT3" 

760  PAUSE 

770  SCALE  0,10,-3,0 

780  M5«i024 

790  ASSIGN*  i  TO  Q5$ 

800  FOR  1=1  TO  99 

810  READ*  1  ;  X<I),Y(I) 

820  PLOT  X(I),LGT(Y<I)) 

830  NEXT  I 

340  PENUP 

850  ASSIGN*  1  TO  » 

860  DISP  "  DO  YOU  WANT  TO  PLOT  A  DATA  FILE    CURVE?  (Y/ 

N>" 

870  INPUT  H% 

880  IF  N$="N"  THEN  GOSUB  970 

890  DISP  "  ENTER  THE  FILE  NAME  OF  THE  DATA  FILE   TO  BE 

PLOTTED.  (NAME;D70i)" 

900  DISP  "  INCLUDE  THE  DISC  DRIVE  WHERE  CAL  FILE  IS  STO 

RED— E.G.  -NAME:D701" 

910  INPUT  Q3$ 

920  M2-1024 

930  GOSUB  1010  !  CALC  THETA  BAR  AND  PLOT  ! 
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940 

IF  P=705 

THEN 

970 

950 

GRAPH 

960 

COPY 

970 

LINETYPE 

i 

980 

CLEAR 

990 

DISP  "PLOTS 

CGrtPLETE 

iOOO  END 

!  "  • 

iOlO  !  READ  AND  PLOT  CAL 

1020  DISP  "ENTER  THE  INDEX  OF  REFRACTION  USED  FOR  THE  C 

AL  RUN.  (AIR=i.O,  WATER=i.3S)" 

1030  INPUT  M 

1040  DISP  "ENTER  THE  D32  TO  BE  USED  IN  COMPUTING  THETA 

BAR  .  " 

1050  INPUT  D 

1060  DISP  "ENTER  THE  DIODE  NUMBERS  FOR  THE  REGION  OF  IN 

TEREST.  (P1,P2)  TYPICALLY  (50,400)" 

1070  INPUT  P1,P2 

1080  H9=D*PI*.025/(L*F*M.) 

1090  P3=H9*(6i .46+Pl) 

1100  P4=H9*<6i.46+P2) 

1110  K6='D*PI/L 

1120  PRINT  -LOWER  LIMIT  OF  THETA  BAR=";P3 

1130  PRINT 

1140  PRINT  "UPPER  LIMIT  OF  THETA  BAR="jP4 

1150  ASSIGNS  i  TO  Q3$ 

1160  FOR  1=1  TO  M2  STEP  4 

1170  READ*  1  }  T(I),I<I) /r(I  +  i),I<I  +  i)  ,T(I  +  2)  ,I(I-»-2)  ,T( 

1+3), 1(1+3) 

1180  NEXT  I 

1190  MOVE  T(1),I(1) 

1200  FOR  I=Pi  TO  P2 

1210  T(I)=K6*T(I) 

1220  IF  I(I)>0  THEN  PLOT  T( I ) ,LGT ( I ( I ) ) 

1230  NEXT  I 

1240  PENUP 

1250  RETURN 

1260  !  LABEL 

1270  Y3=Y/iOOO 

1280  MOVE  10.1,LGT(Y)  »  LORG  2 

1290  LABEL  USING  El$  )  "-",Y3 

1300  MO^>E  -.1,LGT(Y)  9  LORG  3 

1310  LABEL  USING  E2$  ;  Y3,"-" 

1320  RETURN 

1330  FOR  Y=Y  TO  Y2  STEP  E 

1340  E2=LGT(Y) 

1350  XAXIS  LGT(Y) 

1360  IF  RMD(Y,10*E)»0  THEN  1380 

1370  GOSUB  1260  !  LABEL 

1380  NEXT  Y 

1390  RETURN 
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